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Abstract. Given an operad P with a finite Grobner basis of relations, we study the generating functions 
for the dimensions of its graded components P{n). Under moderate assumptions on the relations we 
prove that the exponential generating function for the sequence {dim P(n)} is differential algebraic, and 
in fact algebraic for P is a symmetrization of a non-symmetric operad. If, in addition, the growth of 
the dimensions P{n) is bounded by an exponent of n (or a polynomial of n, in the non-symmetric case) 
I then, moreover, the ordinary generating function for the above sequence {dimP(n)} is rational. We give a 

■ number of examples of calculations and discuss conjectures about the above generating functions for more 

^ ■ general classes of operads. 
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0. Introduction 



We study the generating series for the dimensions of the components of finitely generated operads. 
■ We conjecture (motivated by an analogy with graded algebras, see Subsection I0.2p that, for a wide class 
of generic operads, these series are differential algebraic for symmetric operads and algebraic for non- 
symmetric ones (for symmetric operads, we consider exponential generating functions). We prove that 
this is indeed the case if the operad has a finite Grobner basis and satisfies an additional (mild) condition. 
Moreover, if the dimensions of components of the operad are bounded by an exponential function (or by 
a polynomial function for nonsymmetric case), then the corresponding generating series is rational. We 
also describe several algorithms for calculating these series in various situations, and provide a number 
of examples of calculations. In particular, there are several natural examples of operads for which the 
generating series was not previously known. 



0.1. Main results. The concept of so-called Grobner basis for operads was recently introduced in |DKj . 

l/^ ' A number of important operads admit a finite Grobner basis of the ideal of relations. Here we discuss 

■ the dimensions of components of operads with a finite Grobner basis of relations in the form of their 
] generating functions. 

■ Let V he a finitely generated operad over a field k of characteristic zero. Recall that an exponential 
. . . generating series of V is defined as 

5-1 ■ n^l ' 



We also consider an ordinary generating series 

(0.1.2) Gp(z) := ^ dim 7^(72)2". 



In particular, if P is a symmetrization of a non-symmetric operad CP, then E-p^z) = Gj>{z). 
Our first result deals with non-symmetric operads. 

Theorem 0.1.3 (Theorem 12. 3. ID . The ordinary generating series of a non-symmetric operad with a finite 
Grobner basis is an algebraic function. 

In the more general symmetric case, the analogous result is true under some additional assumptions 
explained below. 

Recall that the leading terms of the elements of a reduced Grobner basis of a symmetric operad are 
shuffle monomials, that is, rooted trees with labeled vertices and leaves. Each shuffle monomial has a 
unique planar realization (see details in Section II. 3p , so we can identify a shuffle monomial with a planar 
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tree whose leaves are enumerated by an initial segment of positive integers permuted by a so-called shujfle 
substitution. We call a set M of shuffle monomials shuffle regular (Definition 13.2. ip if for each shuffle 
monomial m in M and for each shuffle substitution a of its leaves, the monomial m' obtained from m by 
acting with a on its leaves also belongs to M. An operad with a given Grobner basis of relations is called 
shuffle regular if the set of the leading monomials of the elements of the Grobner basis is shuffle regular. 

Recall that a function or a formal power series is called differential algebraic if it satisfies a non-trivial 
algebraic differential equation with polynomial coefficients. 

Theorem 0.1.4 (Corollary I3.2.7|) . Let V be a shuffle regular symmetric operad with a finite Grobner 
basis of relations. Then its generating series E-p is differential algebraic. 

We also consider two special classes of operads which give rise to generating series of more special form. 
The first class consists of monomial shuffle regular operads with an additional symmetry of relations. 
Namely, if the set of relations of a monomial operad V forms a set of all planar realizations of a given 
set of non-planar trees, then the generating series E-p is algebraic (Theorem I3.3.2p . The second class is 
defined by the following bounds for the dimension growth of the components of the operad. 

Theorem 0.1.5 ( Cor ollar ies 1 2 . 4 . 1] and [3l^T|) . LetV be an operad with a finite Grobner basis of relations. 
Suppose that either 

{i) V is non-symmetric and the numbers dimVn are bounded by some polynomial in n 

or 

{ii) V is shuffle regular and the dimensions dimVn o,re bounded by an exponential function a" for 
some a > 1. 

Then the generating series G-p is rational. 
In fact, the growth conditions we need in Theorem 10.1.51 are even a little weaker, see Corollaries 12.4.11 



All our proofs are constructive and provide effective methods of obtaining the corresponding algebraic 
or differential algebraic relation for the generating series. We give a number of methods; depending on 
the situation, one of these may be the most effective. These methods are based on relations between 
right sided ideals in monomial operads, on homological computations, and on using symmetries of the 
relations of operads. For some classes of operads, we derive a general formula for the generating series 
or prove that the generating series is algebraic. We also give a number of examples of calculations. 

Remark 0.1.6. Sometimes there exists an additional integer grading on the operad V such that all 
vector spaces Vn are graded, and such that the grading is additive with respect to the compositions of 
homogeneous elements. Then one can also consider the two-variable generating functions G-p{z,t) and 
Ep{z,t), see Sections \2.1\ and \3.1i All the our above results remain valid after replacing the coefficient 
field Q in the differential and algebraic equations by the ring Q[t]. 

0.2. Motivation. In this paper we study generating series for operads. A lot of examples of such series 
are recently calculated, see [Z]. Moreover, for binary symmetric operads with one or two generators 
(essentially, for quotients of the operad of associative algebras), such generating series has been extensively 
studied for decades under the name of codimension series of varieties of algebras, see \GZ\ IBD] . Here 
we present a new approach to such series based on the theory of Grobner bases for operads developed 



Many results in the theory of operads have analogy in the theory of graded associative algebras, which 
we now partly review. 

For a graded finitely generated associative algebra ^ = © © ^2 © • • • over a field k, its Hilbert 
series is defined as the generating function for the dimensions of the graded homogeneous components: 



Whereas in general the series Ha{z) can be very complicated (e.g., is a transcendental function), for a 
number of important classes of algebras it is a rational function 



and lSXn 



in [DK] . 



Ha(z) = ^(dimk^„)z'^. 



p{z) 



, where p(z), g(z) € 
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Indeed, Ha{z) is rational if A is commutative, or Noetherian with a polynomial identity, or relatively 
free, or Koszul (according to a conjecture of [PPj ). etc. Probably the most general class of algebras with 
rational Hilbert series is the class of algebras with a finite Grobner basis of relations. Finitely presented 
monomial algebras, PBW algebras, and many other natural types of algebras are particular examples of 
this class. For a survey of these and other results on Hilbert series of associative algebras, we refer the 
reader to [U] and references therein. 

We conjecture that the generating series of the most important general examples of operads belong to 
some analogous reasonable class of formal power series that we now describe Let us try to bound this 
hypothetical class. 

Certainly, the generating series of the most useful operads such as operads of commutative. Lie, and 
associative algebras should belong to the hypothetical class of simplest generating series. Since we have 

Ecom = — 1, E^Assoc = ]^ _ ^ 1 E Cie = — hl(l — z), 

we conclude that, at least, the generating series of a generic symmetric operad may not be rational and 
may be exponential or logarithmic. In addition, finite-dimensional operads seem to be simple enough 
to have "general" generating series. Hence the described class should also include the polynomials with 
rational coefficients. Free finitely generated operads are also general; if the generating series of the vector 
space (more precisely, the symmetric module) V of generators of such an operad 3~ is a polynomial 
p{z) = Ey, then the generating series of 3~ is 

where fv{z) = z — p{z) and stands for the inverse function. This means that our class of formal 
power series should also includes algebraic (over Q) functions. 

Note that the sets of quadratic, finitely presented, and other main types of operads are closed under a 
number of operations (such as direct sum with common identity component, coproduct, and composition). 
Thus, it seems reasonable to assume that the set of common generating series should be closed under 
corresponding formal power series operations, that is, the operations that send the pair of formal power 
series f{z) = E-p{z) and g{z) = Eq{z) to 

f{z) +g{z) - z, 

ancfl 

{f-\z) + g-Hz) - z)-\ 

Finally, since the Koszul duality plays an important role in the theory of operads and its applications, 
one would want to have that if a generating series of a Koszul operad T belongs to this class, then the 
series of its Koszul dual J"' should also belong to this class. To ensure this, we assume that the described 
class of formal power series should be closed under the operation which sends f{z) to 

-(r')(-^)- 

Note that while the generating series of Com and Cie are not rational nor algebraic, both these series 
satisfy simple first-order differential equations. In particular, these functions are differential algebraic. 
In view of the consideration above, we state the following claim. 

"'^To prove the equality h{z) — (/^^(z) + g~^{z) — z) ^, where h{z) = -E-p]jq is the generating series of the coproduct 
(=free product with common identity component) of operads V and Q, consider the (minimal) differential graded models 
M-PjMq, and M-pjjg of these three operads. They are DG operads which are free as operads. Let H-p, Hq, H-pY^Q be 
the symmetric modules which minimally generates these three DG operads; then H-pjjQ ~ H-p © Hq. One can consider 
these symmetric modules as complexes (with trivial differential but non-trivial homological grading); let x(^)iX(2)i 
x(^IJ 2) = xi'P) + x(2) be the exponential generating functions of their Euler characteristics. For each operad O, the 
analogous exponential Euler characteristic of its minimal model is 

x{Moy^{z)^ z-x{Ho){z), where x{Mo){z) ^ Eo{z). 

Thus, 

Evuai^) = x(Mpus)"'(^) - z - x{VY[Q) = z ~ x{V) + X(2) = r\z) + g'\z) - z. 
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Expectation 1. The exponential generating series of a generic finitely presented symmetric operad is 
differential algebraic. 

As for the ordinary generating series of non-symmetric operads (i.e., the exponential generating series 
of symmetric operads which are symmetrizations of the non-symmetric ones), the class of such general 
functions should include, at least, the above generating function of Assoc, polynomials over nonnegative 
integers (for finite dimensional operads) and formal power series which are roots of some polynomials 
with integer coefficients (for free operads). In particular, such generating series are algebraic, that is, 
they satisfy some non-trivial algebraic equations over Z[z]. 

Expectation 2. The generating series of a generic finitely presented non-symmetric operads are algebraic 
over Tl\z\ . 

Suppose that the sequence of dimensions dim7^„ of components of an operad has slow growth (e.g., 
dim Vn is bounded by some polynomial of re) . Then in many examples the ordinary generating series G-p 
of the operad V is quite simple. For example, the operad Com has a rational generating series: 

Gcom (z) = Z + + ^3 ^ . . . = 

Such examples suggest the following claim. 

Expectation 3. The ordinary generating series G-p of a generic symmetric or non-symmetric operad V 
is rational if the sequence {dimVn} is bounded by a polynomial in re. 

In particular, if the sequence {dimVn} is bounded by a constant, then it is eventually periodic. 

To support our claims, we consider operads with a finite Grobner bases of relations. The analogy with 
the theory of associative algebras leads us to expect that the generating series of operads with finite 
Grobner bases are in our class. This expectation is suggested by the results of this article. 

0.3. Outline of the paper. In Section [U we give a brief introduction to algebraic operads and Grobner 
bases for them. In particular. Subsections 11.21 and 11.31 contain detailed description of monomials in free 
non-symmetric and symmetric operads respectively. 

In Section [21 we illustrate our ideas in the simpler case of non-symmetric operads. In Subsection 12. 2. if 
we begin to prove Theorem 10 . 1 . 3 1 by finding an algebraic relation for the generating series of an arbitrary 
non-symmetric operad with a finite Grobner basis. Our proof is based on considering of the principal 
ideals with bounded degree (or level) of generators in the finitely presented monomial operad. We provide 
an algorithm to construct polynomial relations of the form yi = F{yi, . . . , y^) for the generating functions 
yi, . . . ,y]y of these ideals. In Subsections 12.2.21 and 12.2.31 we give two other versions of this algorithm 
optimized for the case of few relations and for the case of simple structure of relations, respectively. Since 
the generating series of the operad itself is a linear combination of some yi, we deduce in Subsection 12.31 
that the latter generating series does satisfy an algebraic equation. Then we discuss a bound for the 
degree of the equation. In Subsection 12.41 we investigate the case of operads with subexponential growth 
of the dimensions of its components and prove the first part of Theorem 10.1.51 

In Section [3l we deal with the case of symmetric operads. In Subsection 13.11 we prove a formula for the 
generating series of the shuffle composition of sequences of linear subspaces. This is the first place where 
integral and differential equations appear. In Subsection 13.21 we prove Theorem 10. 1.41 In Subsection 13.31 
we consider a class of shuffle regular operads with additional symmetries among the relations. For these 
operads (called symmetric regular), the sets of leading monomials of a Grobner basis of relations are, 
by definition, closed under arbitrary re-numbering of the leaves. We show that exponential generating 
series of symmetric regular operads are algebraic functions. In Subsection 13.41 we consider shuffle regular 
operads with slow growth of the dimensions of the components and prove the shuffle part of Theorem 10. 1.51 
All provided examples of symmetric operads are collected in Subsection 13.51 

Some remarks and conjectures are collected in the final Section [H First, we give some evidence for 
the conjecture that each binary operad with quadratic Grobner basis has the same generating series as 
some shuffle regular operad (hence, this series is differential algebraic). Second, we discuss some evidence 
showing that operads of associative algebras with polynomial identities are "generic" in the sense of 
our Expectation [1] (and in the sense of Expectation |3l in the case of polynomial growth). Third, we 
remark that the operation of shuffle composition induces a structure of Zinbiel algebra on the ring of 
formal power series. Therefore, it should be interesting to describe the minimal Zinbiel subalgebras 
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which contain generating series of some classes of operads, and to describe the class of "Zinbiel algebraic" 
functions which contains these generating series. 

0.4. Considered examples. We illustrate our methods by a number of examples. For the reader's 
convenience, we give below a brief list of calculations (more or less explicit) of generating series for some 
classes of operads. In most of these examples, the operads are generated by a single binary operation 
(multiplication). As usual, we denote [a,b] := ah — ha and (a,6, c) := {ab)c — a{hc). These examples are: 

• a non-symmetric operad with the identities {{ab)c)d) = and {a{h{{cd)e))) = (Example 12. 2. 

• non-symmetric operad with generalized associativity relation 
xi{. . . ixk-2ixk-i,Xk,Xk+i)) . . . ) = (Example [MJ^; 

• operads of alia algebras, left alia algebras and right alia algebras [Dzhj . that is, the operads 
generated by a single binary nonsymmetric operation defined by the identity {[xi, a;2], X3} -|- 
{[x2,X3],xi} + {[x3,xi],X2} = 0, (where {a,h} := ah + ba), [xi, X2]x3 + [x2, xsjxi + [x3,xi]x2 = 0, 
and xi[x2,X3] +X2[a;3,a;i] -|-X3[xi,X2] = 0, respectively (Example 13. 5. ip . 

• operad J\fV(n of upper triangular matrices of order n over a nonassociative commutative ring 
(Example 13.5.31 for n = 2, Example 13.5.61 for n = 3 and Lemma 13.5.41 for general case). 

• a Lie-admissible operad defined by the set of identities {g{[—, —],..., [— , —])\g S where the 
generator [— , — ] satisfies the Jacobi identity and g is a composition of other generators operation 
$ is a given set of operations on some generators and the additional generator (Example 13. 5. 7p : 

• an operad with the identities [x, [y,z]] + [y, [z,x]] + [z, [x,y]] = and -|- [2:,t][x,y] = 
(Example 13.5. IH as a special case of the previous example). 

Acknov^ledgments. We are grateful to Volodya Dotsenko, Askar Dzhumadil'daev, Samuel Grushevsky 
and Boris Shapiro for stimulating discussions and comments. The main part of the preprint was written 
while the first author was working at the ETH and he wishes to thank all the staff of ETH and especially 
Giovanni Felder who in addition helped to organize a short visit of the second author to ETH. 

1. Operads and trees 

In the brief and rough explanation below, we refer the reader to the books [LV^ IMSS] for the details 
on the basic facts on operads and to [DKj and [LV^ Ch. 8] for the details on shuffle operads and Grobner 
bases in operads. 

1.1. Operads in few words. Roughly speaking operad is a way to define a type, or a variety, of 
(linear) algebraic systems. An operad IP consists of a sequence (or a disjoint union) of sets (vector 
spaces) IP(n) of all multilinear operations (numbered by the amount of inputs) and composition rules 
IP(n) o J'(m) I— 7- y(n + m — 1) which prescribe how one can substitute an input of an operation by the 
result of another one. 

One typically separates two versions of this notion (the so-called symmetric and non-symmetric op- 
erads) depending on whether we allow and do not allow to permute the inputs of operations. The 
permutations of inputs define the action of symmetric groups on operations and compositions, hence, 
the direct computations becomes much harder in the symmetric case. On the other hand, the symmetric 
operads are much more important for applications. The general theory of non-symmetric operads is 
similar but more transparent, since one needs not take care of the symmetric group action. 

Given a discrete set T = U„^iT„ of (abstract) multilinear operations, where T„ consists of the n-ary 
operations, one can define a free operad 3" = 3'(T) = © 3"n generated by T. Here 3"„ is the vector space 

spanned by all n-ary operations which are compositions of the elements of T with each other (and the 
first component 3"i is assumed to be spanned by the identity operation which, formally, does not belong 
to T). Note that in the symmetric case a 'composition' may permute the inputs, so that the set of n-ary 
generators T„ should span a representation of the symmetric group. To separate the non-symmetric and 
symmetric cases, we will denote the corresponding free operads by 3'(T) and J-(T), respectively. 

Each operad can be defined in terms of generators and relations. Therefore, each operad may be 
considered as a quotient of a free operad. That is why we first describe linear bases in the free operads 
(referred to as monomial bases) which are compatible with compositions of operations. The difference be- 
tween the symmetric and the non-symmetric cases is large already on this level. Namely, the definition of 
monomial basis for non-symmetric operads is natural whereas to define monomials for symmetric operads 
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we need an additional structure of shuffle operads introduced in [DK| . We hope that the definitions and 
combinatorics of non-symmetric and shuffle monomials will be clear from the forthcoming subsections. 



1.2. A basis for a free non-symmetric operad: tree monomials. Let us recall the combinatorics 
of trees involved in the description of monomials in operads. 

Basis elements of the free operad are represented by (decorated) trees. By a (rooted) tree we mean 
a non-empty connected oriented graph T of genus zero for which each vertex has at least one incoming 
edge and exactly one outgoing edge. Some edges of a tree have vertex at one end only. Such edges are 
called external. All other edges (having vertices at both sides) are called internal. Each tree has one 
outgoing external edge (= the output or the root) and several ingoing external edges, called leaves or 
inputs. The number of leaves of a tree v is called an arity of a tree and is denoted by ar(v). 

Consider a set of generators T = Un^iT^ of a free operad. We mark a tree with a single vertex and n 
incoming edges by an arbitrary element of T„ and we call such tree a corolla. We say that a rooted tree 
T is internally labeled by the set T if each vertex is labeled by an element from T such that the number 
of inputs in this vertex coincide with the arity of corresponding operation. 

A rooted tree is called planar if for each vertex there is an ordering of all incoming edges. (Therefore, 
there is a canonical way to project a tree on a plane.) 

Proposition— Definition 1.2.1. A canonical basis of a free non-symmetric operad 3"(T) is enumerated 
by the set B(T) of all planar rooted trees internally labeled by the elements of the set T. We shall refer 
to the elements of this basis as planar tree monomials. 

There are two standard ways to think of elements of an operad in 
terms of its generators. The first way in terms of tree monomials rep- Figure 1. A pla- 

resented by planar trees and the second one is in terms of compositions nar tree monomial 

of operations presented by formulas with brackets. Our approach is 
somewhere in the middle: in most cases, we prefer (and strongly en- 
courage the reader) to think of tree monomials, but to write formulas 
required for definitions and proofs in the language of operations since 
it makes things more compact. 

A particular example of a planar tree monomial is presented in 
Figure [TJ The corresponding operation may be written in a following 
way g{f {/{-,-), g{-, -,-)),-, f {-,-)). It is also convenient to use vari- 
ables/letters for the inputs of operation. In a non-symmetric operad, 
we enumerate the inputs from the left to the right, so that the input variables are (from the left to the 
right) xi,X2,X3, . . . For example, the operation from Figure [1] will be written in the following form: 

9{fif{xi,X2),g{xs,X4,X5)),Xe,f{x7,X8)). 

The composition rules for tree-monomials are given by the concatenation of trees (or composition of 
operations). The notion of divisibility is defined in the following way: A monomial v is divisible by w if 
there exists a planar subtree of v isomorphic to tt; as a planar labeled tree. (Notice that the root of w 
not need to coincide with the root of v.) For example, the tree-monomial that represents the operation 
f{-,g{-,-,-)) is a divisor of a tree-monomial drawn in Figure [T] which corresponds to a unique minimal 
subtree which contains the internal edge going from the ternary operation g in the upper level to the 
binary operation / in the intermediate level. 

Let us specialize two cases of divisibility. A monomial w of some arity m is a left divisor of v if 
V = w{vi, . . . , Vm) for some monomials vi, . . . , Vm- By the other words, ri; is a left divisor of v if there is a 
planar subtree w' in v which is isomorphic to vo and has the same root as v. Similarly, a monomial u of 
some arity / is a right divisor of v if f = vo{xi^ . . . , u(xi, . . . , Xj_|_;_i), rEj_|_;, . . . , x;+m-i) (where u is in ith 
place) for some monomial w of some arity m. By the other words, u is a right divisor of v if all leaves 
of the subtree u' in v isomorphic to u are also leaves of v. For example, the empty monomial and the 
monomial v itself are both left and right divisors of v. 

1.3. A basis in free symmetric operad: shuffle monomials. The free symmetric operad may not 
have a monomial basis compatible with all possible compositions. The way to avoid this problem is to 
construct a basis compatible with some compositions. A class of compositions we would like to preserve 
is called shuffle compositions. A collection of all multilinear operations with the prescribed rules for 
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shuffle compositions form a shuffle operad ( |DKj ). Any symmetric operad V may be considered as a 
shuffle operad denoted by . Reversely, one can recover from all information concerning V but the 
action of the symmetric groups. In particular, the generating series of the operad V and are the same. 
All definitions and conclusions below are given for shuffle operads and their monomial bases. However, 
symmetric operads are close to the nature, hence all examples of shuffle operads discussed below are 
examples of symmetric operads considered as shuffle ones. 

Let us define a basis in the free symmetric operad J^(T) compatible with the shuffle compositions. The 
elements of the basis will be called shuffle monomials. The main difference from a non-symmetric world 
is that one should have an external labeling of an operation (i.e., inputs should be ordered). We say that 
a rooted tree with n incoming edges (=leaves) has an external labeling if the set of leaves is numbered 
by distinct natural numbers from 1 to n. 

Proposition— Definition 1.3.1. The basis of a free shuffle operad J~{T) generated by a set T is numbered 
by the set B{T) of all rooted (nonplanar) trees internally labeled by the set T and having external labeling. 
We shall refer to the elements of this basis as shuffle tree monomials or simply shuffle monomials. 



Note that whereas the shuffle tree monomials are originally considered as abstract graphs, we would 
need their particular planar representatives that we now describe. 

In general, an embedding of a (rooted) tree in the plane is determined by an ordering of inputs for 
each vertex. To compare two inputs of a vertex v, we find for each of the inputs the minimal leaf that 
one can reach from v via it. The input for which the minimal leaf is smaller is considered to be less than 
the other one. The corresponding embedding of a shuffle monomial to the plane is called its canonical 
realization (in JDK', 3.1], it is referred to as canonical planar representative) . Thus we identify a shuffle 
monomial with its canonical realization. Note that a canonical realization is a planar tree whose leaves 
are enumerated by an initial segment of positive integers. The leaves may be permuted by substitutions 
of a particular kind; we call them shuffle substitutions. The only assumption is that it should preserve 
the ordering of minima of the leaves in each vertex. 

Similarly to the non-symmetric case there are two languages to think about the tree-monomials and 
operations in the free operad. We still recommend the reader to think on shuffle tree-monomials as 
labeled planar trees but most of the formulas in examples are written in the more compact language 
of operations. See Figure [2] below for the comparison of this two languages, with the monomials 
9{f{fixi,X3),g{x2,fix4,X9),g{x5,XQ,xu))),X7,f{x8,xio)) and /(xi, 5(x2, X3, X4)) represented as trees. 

Let us explain the notion of divisibility of shuffle monomials. We say that a monomial v is divisible by 
a monomial w if there exists a labeled subtree u inside v isomorphic to w as labeled tree and such that 
the local order of the minima of the leaves that can be reached to the inputs of a subtree u inside v is 
isomorphic to the order of leaves in w. In other words, if we attach to the input labeled hy i in w the 
minimum of the leaves in the corresponding subtree of v (the subtree coming to this input in u C v) 
the following inequalities should be satisfied: i < j 



Hi < Uj. 



an 



example of divisibility of shuffle monomials. 
Figure 2. Divisibility of shuffle monomials 



Figure [2] provides 
In the shuffle mono- 
mial represented by 
the left tree we 
encircle a subtree 
with one internal 
edge. This sub- 
tree is a divisor of 
the left monomial. 
The tree in the 
middle is the encir- 
cled divisor where 
the leaves are num- 
bered by the min- 
ima of the corre- 
sponding subtrees. 

The right tree represents a shuffle monomial where we put the subsequent numbers on leaves according 
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Let us define two particular cases of divisibility. We say that a shuffle monomial v is right divisible 
by a shuffle monomial w if there v is divisible by w and, in the notation above, all external vertices of 
the subtree u are also external as vertices of v. This means that there exists another shuffle monomial t 
such that the monomial v is obtained from t by replacing one of its input (say, the m-th) by a subtree 
isomorphic to u. In this situation, we call the monomial v a shuffle composition of t and w and write 
V = to^ w. Note that this notation is insufficient to define the composition uniquely because the former 
should contain also an information on "shuffling" substitution of the inputs of t and w [DK^ Prop. 2]. 

Analogously, a shuffle monomial u (say, of arity m) is a left divisor of v if there is a labeled subtree u 
inside v with common root with v which is isomorphic to w as a labeled tree and with the local ordering 
of the inputs with the same property as above. In this case one can represent u as a multiple shuffle 
composition 



for some shuffle monomials fi, . . . ,Um- We denote this multiple composition simply by u{vi, . . . ,Vm)sh- 
Again, this notation is insufficient since it does not contain a part of information on the permutation of 
the inputs of v. 

1.4. Monomial operads and Grobner bases. Let T be a set of generators of a non-symmetric or 
shuffle operad V, that is, each element of "P is a linear combination of compositions of elements of T 
with each other and the identity operator. Each component Vn has a natural structure of the quotient of 
the component J-"(T)„ of the (non-symmetric or symmetric) free operad J~{T) by some vector space Cn- 
Then the suboperad C = U„j.iC„ in J^(T), being a kernel of the natural surjection J-{T) — )• V, forms an 
ideal in J^(T) called the ideal of relations of V. The operad V is called monomial if each vector space C„ 
is spanned by monomials (non-symmetric or shuffle, according to the type of the operad V). In this case, 
the ideal C is called monomial as well. For a monomial operad V, all monomials which do not belong to 
Cn form a linear basis of Vn for each n. 

Given an order on monomials compatible with compositions (for the discussion on admissible orderings, 
see [DKl §3.2]), one can define a leading monomial for each element of the free operad. For the ideal C of 
the relations of an operad V as above, let C be the span of the leading monomials of all its elements. Then 
C is a monomial ideal such that the monomial basis of the quotient operad V = J-{T)/C is also a basis of 
V (if one identifies the monomials in -F(T) with their images in V). For instance, the dimensions of the 
corresponding components of V and V are the same, so that these two operads have the same generating 
series. A Grobner basis of the ideal C is a set of elements of C such that their leading monomials generate 
the monomial ideal C |DKl §3.5]. Grobner bases gives a way to construct the monomial operad V starting 
from V. This leads to a linear basis of V and to the dimensions of its components. 

Note that if an operad V admits a finite Grobner basis then the generating series of V is the same 
as the generating series of the operad V defined by a finite set of monomial relations (i.e., the leading 
terms of the elements of the Grobner basis). In this paper, we study the generating series of such operads 
v. The operads which admits finite Grobner basis includes PBW operads ([H], |DKt Cor. 3]), operads 
coming from commutative algebras |DKt §4.2] and many others. Some new nonquadratic examples are 
presented below. 

We split the main part of this article into two sections. Section [2] contains the results for generating 
series of the non-symmetric operads. Section [3] contains a (somewhat more complicated) version for 
symmetric operads. In the additional Section |4] we discuss some conjectures and connections with other 
types of operads. 





2. NONSYMMETRIC OPERADS 



2.1. Generating series and compositions. Suppose that the subset M C B(T) defines a monomial 
basis of a non-symmetric operad 7 := 3~(T)/(<I>) (where ($) denotes the ideal in the free non-symmetric 
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operad 3"(T) generated by the set $ C 3"(T) called the set of relations of CP). The (ordinary) generating 
series of CP is defined as the generating function of the dimensions of its components as follows 

Gt{z) := ^ dim CP(n)z" = ^ z^'^^l 

Moreover, if there exists an additional grading of the set of generators T such that all relations from 
$ are homogeneous with respect to this internal grading, one can consider in addition a generating series 
with two parameters 

Gj>{z,t) := ^dimiCP(n)z" = ^ tl^lz'^K^) 

On the other hand, we will omit the parameter t if the internal grading is defined by the arity of operations. 
For example, if all generating operations has degree 1 with respect to internal grading and has same arity 
fc > 1, then all operations of arity n(k — 1) + 1 should have the same grading n. 

Let us consider in details the case of one composition. Namely, let /i be an m-ary generator in the 
free operad 9'(T) (i.e., fj, S T^), and let CPi, . . . ,Tm be some graded vector subspaces in 3"(T). A non- 
symmetric composition /x(CPi, . . . , Tm) of vector spaces Ti, . . . , CP^ is the subspace of 3"(T) whose n-th 
component is spanned by all possible compositions 

fi{pi, . . . ,pm), where pi G Tiih) and '^ki = n. 

Lemma 2.1.1. The generating series of a non- symmetric composition of vector spaces /x(CPi, . . . , CP^) is 
the product of generating series of inputs: 

(^, t) = t'^' (Gn iz,t)-...- Gy^ (z, t)) . 

Proof. It is obvious from the description of the canonical basis in a free non-symmetric operad that 

m as a graded vector space. D 

2.2. System of equations for generating series. We give here a result which is a key point of 
Theorem 10.1.31 It gives a system of algebraic equations which will lead (via elimination of variables) to 
the algebraic equation for Gy, see Subsection 12.31 The proof of this result given here contains some core 
algorithms enabling computations in particular examples. 

Theorem 2.2.1. For a given non-symmetric operad CP with a finite set of generators and a finite Grobner 
basis there exist an integer N and a system of algebraic equations on N + 1 functions yo = yo{z,t),. . . , 
VN = yN{z,t) 

(2.2.2) y,=t^^ Yl li-Vs.- ■■■■Vs,^ fori = l,...,N, 

s&[Q..N]''i 

such that Gy{z, t) = YliLo Uiiz, t), yo = z and yi{0, t) = yi{z, 0) = 0, i = 1 . . . N . The numbers q\ G {0, 1} 
and the nonnegative integer numbers di, and N are bounded from above by some functions of the degrees 
and the numbers of generators and relations of the operad CP. 

Note that, under the initial conditions the solution of system (j2.2.2p is unique. It follows that one can 
consider ()2.2.2p as a system of recursive equations on the coefficients of series yi . We will refer to ()2.2.2p 
as the system of recursive equations. 

Below we present two different proofs of this theorem (in Subsections 12. 2.1] and [?.2.2I below). It is 
better to say that we present two different algorithms to construct the above system of equations. In 
addition, we give also an idea of a third algorithm. 

The first one is more efficient in the case of many relations of relatively low arity. Additionally, the 
first proof gives the positivity of coefficients ql which we use in Corollary 12.4.11 The second algorithm 
seems much more useful in case of few relations. Namely, the number of additional variables is 
typically much lower than in the first algorithm. On the other hand, the second algorithm do not 
allow to have only positive coefficients ql, that is, it gives a slightly weaker version of Theorem 12.2.11 
Moreover, in Subsection 12.2.31 we present an example of computations based on a third idea, that is, we 
use combinatorics of homology of the operad. We hope that a reader can get the general idea of this 
method from the given example in order to compute the functional inverse to the generating series of a 
given operad. 
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Note that the systems of equations similar to the system 12.2.21 above are sometimes appear in combi- 
natorics and computer science in symbolic methods for combinatorial structures (see e.g. |FS] ) as well as 
for the context-free languages (see e.g. |ChSch| ) . 



2.2.1. First proof of Theorem V2.2.1\ 

Proof. Suppose that an operad 7 has a finite set of generators T and a finite set of monomial relations $. 
(It suffices to consider monomial relations since we are dealing with generating series. Therefore, there 
is no difference between the relations that form a Grobner basis and corresponding monomial relations 
presented by the leading terms of the Grobner basis.) Let d be the maximum level of leaves of elements 
of $ (by the level of a vertex/leaf in a tree we mean the number of the edges in the path from the root 
to this vertex/leaf). As mentioned in Section 11.2.11 everv monomial v in the free operad 3"(T) can be 
identified with a rooted planar tree whose vertices are marked by elements of T. Given such a monomial 
V, by its stump b{v) we mean its maximal monomial left divisor such that the leaves of b{v) have levels 
strictly less than d. In other words, b{v) is the submonomial (rooted subtree) of v which consists of the 
root and all vertices and leaves of v of level less or equal to {d — 1) and all edges connected them. 

Let Stump be the set of all stumps of all nonzero monomials in !P. Let N be the cardinality of this set. 
The elements 6i, ... ,67V of Stump are partially ordered by the following relation: 

bi < bj iS i ^ j and bi is a left divisor of bj . 

Let Mfe. be the set of all monomials in (= the monomial basis of) the right-sided ideal (bi) C 7 generated 
by bi, and let 

j:bi<bj 

Then the pairwise intersections of the sets M;,. are empty. Moreover, the disjoint union IJiLi ^bi is the 
monomial basis of the operad 7. We have 

N 

Gt{z) = ^yi(2;), 
1=1 

where yi{z) = an(M')^'^^ generating series of the span of the set Mb^. For every element 

(=operation) /i G T of some arity n, let us define the numbers jf_i{ii, . . . , in) for all 1 ^ zi, . . . , ^ as 
follows: 

f 0, if u(M^,...,M^) = Oiny 
(2.2.3) jf,{ii,---,in) = { ' ~ - — ' 

[ j, if the stump b{fi{Mb^_^ Mb^J) = bj. 

Note that the nonzero sets of the type /u(Mfe.^, . . . ,^6,^) have empty pairwise intersections. Let v 

be a nonzero monomial in CP with the root vertex labeled by fi. Then v G fi{M,b-^, . . . ,Mf,.^) for some 

fi,ii, . . . ,in, that is, v = fi{vi-^, . . . ,Vi^) where the monomial subtrees Vj, G Mf,. . are uniquely determined 

by V. Hence, v £ "Mb^ where j = jf_i{ii, . . . ,in) from (I2.2.3p . As soon as the degrees of the relations are 
less than or equal to d we come up with the following disjoint union decomposition for all j = 1 . . . A^ 

Mb- = ly ^(Mfc.^ , • • • , ^bi^ ); where fj, is the root vertex of any v G Mb- . 

j^{ii,...,i„)=j 

This equality implies equations ()2.2.2p for the generating functions yi{z) = G^ an(M')^'^^' ^ 

Notice that if a stump bj do not have leaves of level {d — 1) then the corresponding set Mb- consists of 
one element bj. For example. Mi always consists of one element representing the identity operation and 
Gmi = z. In fact, the substitution i/q = Gjyti = z is already made in the system (I2.2.2|) . Therefore, one 
can reduce the number of algebraic equations in the system (I2.2.2p to the number of the stumps of level 
{d — 1) which contain at least one leave of level (d — 1). In particular, if the operad is a PBW operad 
then the number of recursive equations is one greater than the number of generators. 
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Example 2.2.4. Let Assoc be the operad of associative algebras considered as a non-symmetric op- 
erad. Namely, Assoc is generated by one binary operation subject to one quadratic relation 

fi{fi{a,b),c) = fi{a, fi{b,c)). It is well known that this relation forms a Grobner basis according to 
the standard lexicographical ordering of monomials \DK.\ [H]. The set Stump of stumps of level less or 
equal to 1 consists of the identity operator 1 and the unary operation fi. So, there are two sets of the 
type M, that is. Mi = {1} and M^. Thus, we get the following system of equations: 



Gj^ssoc — z -\- y^, 
yii = z"^ + zy^ 



G^ssoc — Z -\- zG J^ssoc G j\^sgQc(yZ) 



2.2.2. Decreasing the number of equations. Now we present one more algorithm that allows to derive the 
system of equations (|2.2.2p . The advantage here is that in many examples the number of equations is 
much less than in the algorithm given in Subsection 12.2.11 above. On the other hand, the coefficients q\ 
might be negative (i. e., q\ G {—1,0,1}). One may consider this algorithm as a generalization of the 
algorithm for binary trees with one relation presented in [Rj. 

The key point of the algorithm is the following. In the notation of Subsection 12.2. It this new algorithm 
gives a system in the generating series of the sets M^. instead of M;,. . Then Gp is equal to the sum of 
z and those Gm;,. where bi runs the generators of the operad T. Therefore, to find Gp it is sometimes 
sufficient to solve only a part of the system which allows to express these new variables. Thus, the 
number of equations (=the number of variables) of the reduced system can be less than the number N 
of equations (and variables) of the system constructed in Subsection I2.2.1[ 

Second proof of Theorem \2.2.1\ Let J" be as above, i.e. IP is a finitely presented (non-symmetric) mono- 
mial operad with a finite set of generators T and a finite set of relations We suppose that the set of 
relations is reduced, namely, g is not divisible by g' for every distinct pair g^g' G Let M be the set 
of monomials in the free operad 3'(T) which are not divisible by the relations from In other words, 
M is a monomial basis of J". We also suppose that the unary identity operation 1 belongs to the set of 
monomials M. We call it the trivial monomial. 

Consider a free operad 3'(T) and the set of all tree monomials B(T) generated by the same set T. With 
each monomial v S B(T) one can associate the subset 3"^ C B(T) consisting of those basis elements that 
are left divisible by t;. In other words, the set 3"^ is a monomial basis of a right ideal f o 3" generated by 
V. Obviously, for the identity operation 1, we have 3"i = B(T). Each collection {vi, . . . ,vi} of monomials 
defines the left common multiple denoted by U . . . U vi]. The monomial [viU . . . U vi] is defined as 
the smallest element in 3"^^ n . . . n 3"^,; . The left common multiple (if exists) should be the unique tree 
given as union of its subtrees Vi. As the definition of the left common multiple one may use the following 
identity: 

3"^i n . . . n 3"^; = 3"[„j^u...ut)i] 

For the case of empty intersection 3"^,^ n . . . n 3"„, we set [viU . . .U vi] to be zero. 

With each monomial f G M we associate the monomial basis M^, of the corresponding right ideal voM,. 
Obviously, we have = M H 3"^, . 

Consider a given nontrivial monomial v G M. Suppose that the root generator of is a /c-ary operation 
/X, that is V = fio [v^^ . . . , v'') where denotes the subtree which grows from the i-th incoming arrow of 
the root vertex of v. Let be the subset of the set of generating relations $ such that the corresponding 
relations have a nontrivial left common multiple with v in the free operad 3"(T), i.e. 

g e<i>v C <^ <^[vU g]^0 &i g e <^> ^%nJg^O k g 

If the left common multiple [v U g] is non-zero then the fi should be a left divisor of v, hence, we have a 
decomposition g = fio [g^ ^ . . . ^g'^) for some g^,. . . ^g'^ . We get the following recursive relation for the set 
M^: 



(2.2.5) =//o(M„i,...,M^fe)\ 



Identity (I2.2.5P should be clear from the observation that each monomial from ^ o (M^i , . . . , M^^ ) either 
belongs to M^, or is divisible from the left by the relo-tion q G 

11 



Denote by yv{z) (or yv{z, t) if the operad is Z^-graded) the generating series of the span of the set My. 
Combining the formula ()2.2.5p with the inclusion-exclusion principle, we count the cardinalities and the 
generating series of the sets M^, and their intersections. We have 

/ \ 



(2.2.6) 



t 



E 

s=0 



E 



{9l,--;as}C<S>v- 

\[vUgiU...Ugs]j^O 



y[viUglu...Ugl] ' • • • ' yKUgJ-'U-UgJ] 



where (as above) fi is the generator placed in the root vertex of both v and gj , and the subtrees growing 
from the i-th incoming arrow of fi are denoted by f * and g^ respectively. The internal summation is taken 
over the subsets {gi, . . . , gs} C of cardinality s such that the corresponding left common multiple 
[v L) gi L) . . . L) gs] is different from zero and is not divisible by any element of <I>. In order to finish the 
proof of Lemma 12.2.11 we have to explain the finiteness of the number of the recursive equations (j2.2.6p . 
This is made in the next Lemma. 

Lemma 2.2.7. There exists a minimal finite set of monomials r(J') C M satisfying the following con- 
ditions. 

• The identity operator 1 and all generators T belongs to T(CP). 

• Suppose that v = fJ.{v^, . . . ,v^) £ T(y) and a collection {gi, . . . ,gs} C has a nontrivial left 
common multiple w = [v L) gi L) . . . L) gs] ^ with v (in particular, each gj has the form gj = 
fJ-igj, ■ ■ ■ , gj) )■ Then for each i = 1, . . . ,k the monomial [v^ U glL) . . .U gl] either belongs to T(J') 
or is divisible by a relation from 

Proof. Suppose that d > 1 is the maximum of the levels of the relations in $. Since the level of the left 
common multiple [v U w] is bounded from above by the maximum of the levels of v and w, we conclude 
that the level of any monomial from T(CP) should be strictly less than d. The set of monomials with the 
bounded level is finite thus T(CP) is finite. □ 

In fact, this bound is quite large. We will see in examples below that the cardinality of T is much less 
than the number of monomials of height not greater than {d — 1). 
It is clear that 



(2.2. 



yi 



At6T 



Therefore, for v G T(CP) we get a system of algebraic equations (|2.2.8p and ()2.2.6p using a finite set of 
unknown functions {yv\v £ □ 

Example 2.2.9. Let Q be a non-symmetric operad generated by one binary operation (, ) satisfying the 
following two monomial relation of arities 4 and 5: 

(2.2.10) ((a, b),c),d) = and (a, (6, ((c, d), e))) = 0. 

The set T(Q) from Lemma 12.2.71 consists of the following 5 elements: 

r(Q) := {1; (a6); ((a6)c); {a{{bc)d)); ((a6)((cd)e))}. 
The corresponding system of recursive equations is 



(2.2.11) 



y{ab) ) 



yl 



y{ab)yi 



yi = z 

y(ab) = 

y{(ah)c) = yiy(ab) 

y(a((bc)d)) 



yiy{{ab)c) + y{(ab)c)y{a{(bc)d)): 

y(ab)y{a({bc)d)) + y{(ab)c)y{a{{bc)d))-, 



y{{ab)c)yi 

yiy({ab)c) -y{{ab)c) - yiy{{ab){{cd)e)) +y{(ab)c)y{iab){(cd)e)), 
y({ab){{cd)e)) = y(ab)y{{ab)c) - yf(ab)c) - y(ab)y{(ab){(cd)e)) +y({ab)c)y{(ab){{cd)e))- 



To exclude additional variables we make the following linear change of variables: y := yi, fs := yi—y({ab)c)^ 
va ■■= yi - y(a{ibc)d)), V5 ■= y{{ab)c) - y({ab)({cd)e))- The systcm (|2.2.1ip is then equivalent to the fohowing 

y - z = V3V4,, 

y-V3 = {V3- Z)V4 =f- V3 - Z = ZV4, 

y-vi = V3V5 

y-V3-V5 = {V3- Z)V5 V3 - V4 + V5 = ZV5. 
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(2.2.12) 



y - Z = z{l + V4,)V4, 
y -y^ = Z{1+ Vi){Vi 

Gq{z) = y = z 



(1-2-^2)2 



2.2.3. Computations via homology. In this section we show how one can simplify computations in some 
cases using the monomial resolutions of operads with finite Grobner bases introduced in [DKl]. First, 
we remind a description of a monomial basis in this resolutions. Second, we present a particular example 
where such a description allows to compute the generating series. The corresponding computation using 
the first two methods (given in two previous examples) became extremely hard compare to what the 
homological method can give. Theoretically, the monomial description of a resolution from [ DKlj allows 
to produce an algorithm similar to the one given in Section 12.2.11 starting from a given finite Grobner 
basis. But the involved combinatorics became extremely hard as soon as the complexity of intersections 
of leading terms of monomials grows up. Therefore we decided not to give all the details of this algo- 
rithm. However, a couple of examples given for the non-symmetric (Example I2.2.14p and the symmetric 
(Example I3.5.7P cases should convince the readers that in real computations they should not close their 
eyes on this method. 

Suppose that i? is a free resolution (=DG model) of an operad T generated by a differential graded 
vector space Q. (In other words, R ~ 3'(Q) as an operad and there exists a differential d on R such that 
the homology operad of R is isomorphic to J") . Then the generating series of J" and Q are related by the 
equality 

(2.2.13) Gyiz) - GQiGj.iz)) = Gy{z) - Gy{GQiz)) = z, 

where Gq{z) is the generating series of the Euler characteristics x of the components Q: 

GQ(z)=^x(Q(n))z". 

Let CP be a finitely presented operad with a set of generators T and a set of monomial relations Let 
us recall a basis in a free monomial resolution of the operad CP. 

quasi 

Proposition (see [DKl] ). There exists a free resolution {R,d) CP such that the set of free generators 
of R consists of the union of the set T and the set CK elements of which are numbered by the following 
pairs: a monomial v £ B(T) and a set {wi, . . . ,Wn} of labeled subtrees (=submonomials) of v satisfying 
the following two conditions. 

[hi) Each Wi is isomorphic to one of the elements of ^ as a planar labeled tree. 

(/i2) Each internal edge of the monomial v should be covered by at least one of subtrees Wi. In other 
words, there is no decomposition v = v' ov" such that each Wi is a subtree of v' or a subtree of v" . 

The homological degree of each generator x £ T is set to be zero and the homological degree of the 
generator {v, {wi , . . . , Wn}) G CK zs set to be k. 

Note that the monomial f in a pair {v, {wi , . . . , Wn}) € ?C is uniquely defined by the set of submonomials 
{wi, . . . , Wn}- For such a pair we will use a notation W. 

Let us present an example where we use this description of a basis in a monomial resolution to get a 
functional equation for the generating series. 

Example 2.2.14. Consider a non-symmetric operad generated by one binary operation (,) which 
satisfies the following relation of degree A; > 2 (i. e., of arity A; + 1): 

rfc := (xi, (. . . , (xfc_2, (xfc_i, Xfc, Xfc+i)) . . . ) = 0, 

where (o, 6, c) denotes the associator (a6)c — a(6c). 

The operad Q2 defined in the same way is Assoc. So, one might consider the identity = of for 
A; > 2 as a weak version of associativity. 



Proposition 2.2.15. The unique generator ri^ of the ideal of relations forms a Grobner basis of relations 
in Qfc. 

13 



Proof. Let us show that all possible s-polynomials reduces to zero. Indeed, there are exactly {k — 1) 
intersections of the leading terms of the relation (r/j): 



k-2 1-2 k+l-l 

Sk+l+2 ■■= {xi{x2 (. . . {{xk-iXk){xk+i (. . . ixk+i-ii{xk+iXk+i+i)xk+i+2 )) • • •) for / = 1, . . . , - 2 

and 

k-2 fe-4 fe+1 k-2 

S2k ■= {Xl{x2{- ■ ■ (. ..{x2k-4{{x2k-3X2k-2)x2k-l)-- ■) X2k) ■ ■ ■) ■ 

The lower index in the left hand side corresponds to the number of leaves/inputs in a monomial. The 
corresponding s-polynomials looks as follows: 

Sk+l+2 {Xl{. . . {{Xk+lXk+l+l)Xk+l+2) •••))- ■ ■ {{Xk-lXk)Xk+l{. . . {Xk+l+lXk+1+2) •••))•• •)) 

S2k ^ • • {Xk-l{{Xk{- ■ ■ {x2k-A{{x2k-3X2k-2)x2k-l)) ■ ■ ■))x2k)) ■ ■ •)) 

- (Xi(. . . {Xk-2{{Xk-l{- ■ ■ {x2k-3{x2k-2X2k-l)) ■ ■ ■))x2k)) ■ ■ ■))■ 

It is easy to see that each monomial of the form 

. . {Xk-2{f{Xk, ■ ■ ■ ,Xj)) ...)), 

where j ^ k + 1 and / is an arbitrary iterated composition of the operation (, ), is reducible via Vk to the 
monomial 

rrij = . . {xj-2{xj-i,Xj)) ...)). 

Thus, both monomials in each s-polynomial above are reduced to the same monomial ruj for suitable 
choice of j. This means that each s-polynomial is reduced to ruj — rrij = 0. □ 

As soon as the Grobner basis is chosen it remains to compute the generating series of homology for 
the corresponding monomial replacement. 

The operad Qfc is generated by one binary operation. Therefore, the tree-monomials under consider- 
ation are rooted planar binary trees where all internal vertices are labeled by the same operation (,). 
Therefore, we will omit this labeling with no loss. In order to specify the tree-type of monomials we say 
that the monomial operation {{xiX2)x^) corresponds to a planar rooted binary tree with two internal 
vertices and one internal edge that goes to the left and the monomial operation (a;i(a;2a;3)) corresponds 
to a binary planar rooted tree with one internal edge that goes to the right. In particular, the leading 
term of the unique element {vk) of a Grobner bases corresponds to a planar rooted binary tree of level k 
with k internal vertices and the unique path that contains all {k — 1) internal edges of a tree. This path 
starts at the root vertex and first have {k — 2) edges that goes to the right and the last edge goes to the 
left. 

For any given element (f , {u)i, . . . , G % all wi should be isomorphic to the leading term of the 
relation Vk as a planar binary tree. Therefore, there exists exactly one submonomial which contains 
the root vertex of v. Without loss of generality we assume that this submonomial is wi. The pair 
(f , {w2, ■ ■ ■ iWn}) will no longer satisfy the property {h2) of the elements in 5{ but can be presented as 
the composition of a generator (, ) taken several times and pairs from CK. In other words, one has to 
present the decomposition of the set of submonomials {1112^ ■ ■ ■ ,Wn} into a disjoint union of subsets such 
that Wi and Wj belong to the same subset if and only if there exists a submonomial v' in v and a subset 
{wi, tUj, . . .} C {w2, . . . , Wn} such that a pair {v' , {tUj, Wj,. . .}) is isomorphic to a pair from 3i. Let us 
show that in the case of operad this decomposition contains at most two subsets. Indeed the unique 
left internal edge of a submonomial wi may not belong to any other submonomial Wi for i > 1. Consider 
a decomposition of a monomial v = v' o v" according to this left internal edge. Then each Wi should be 
a submonomial of v' or v" and we have a decomposition 

(2.2.16) {W2, . . .,Wn} = {W„^2),-- ■,Wa(^l)} U . . .,W„(^ri)} 

for an appropriate permutation a. There also might be the cases when the first or the second subset is 
empty. There exists a pair of uniquely defined (probably empty) submonomials vi C v' and V2 C v" such 
that both pairs (vi, {u;ct(2), • • • j i''^2,{'Wa{i+i), ■ ■ ■ ,u)a{n)}) isomorphic to the elements of 

!K. Moreover, the monomial tree v can be uniquely presented as a composition of several generators and 
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monomials vi and V2- We get a recursive formula for the generating series y of the Euler characteristics 
of elements in !K: 

1=2 1=1 

Here the first summand corresponds to the empty set {if 2; • • • ; Wn}; the second summand corresponds to 
the empty set {^^(i+i), . . . , u'o-(n)}; the index / corresponds to the number of internal vertices in a subtree 
v' which do not belong to none of the vertices of submonomials in {W(j(2) ) • • • i 'Wa{i)}'j the third summand 
deals with both nonempty sets in decomposition (j2.2.16p ; and the fourth summand corresponds to the 
empty set {u;o-{2) , • • • , 'U^(7(/)}- The powers of z are equal to the number of leaves coming from the internal 
vertices that do not belong to submonomials vi and V2- The minus sign comes from the homological 
degree since we remove exactly one element wi from the set of submonomials. 

Finally we have the following quadratic equation for the functional inverse series Gqj.(z)^^ = z — z^—y: 

which is equivalent to the algebraic equation of degree {k + 2) on the generating series Gq^. . 

2.3. Single algebraic equation for generating series. The classical elimination theory implies ex- 
istence of algebraic equation on a function Gy from the system (I2.2.2p (see explanation below). See 
also |ChSch] and the appendix B.l in [FS| and references therein where the same theorem is proven for 
context-free specifications and languages. 

Theorem 2.3.1. The generating series Gj> of a non- symmetric operad CP with a finite Grobner basis is 
an algebraic function. 

Starting from the system (I2.2.2p it remains to use the following 

Lemma 2.3.2. Suppose that the formal power series fi, ■ ■ ■ , fn without constant terms in variables t and 
z over Q satisfy a system of algebraic equations of the form 

fif"' =9i{fl,---Jn) 

for each i = 1,. . . ,n, where gi is a homogeneous polynomial in n variables of degree di ^ 2. Then the 
power series fi satisfy a polynomial equation 

Qif) = 0, 

where Q is a non-constant polynomial with coefficients in Q[t,z] such that degQ ^ (nr=i ^«)^- 
Proof of Lemma \2. 3. 2\ The above system has the form 

F = G{P), 

where F = (/i,... and G = {t'l'^gi, . . . ,t-^'"gn), or 

H{F) = 

with H = ld —G. Note that the Jacobi matrix J = dH/dF is non-degenerate, because det J = 1+0{F) ^ 
0. Let K be the field of rational fractions Q(t, and let L be some its algebraically closed extension 
which contains the ring of formal power series in t and z. Obviously, the variety ^ C of the solutions 
of the above system is 0-dimensional. 

Therefore, there exists a non-trivial polynomial T{x) over L such that T{fi) = 0. By Bezout's theorem, 
one can take T such that degT ^ 0^=1 effective Hilbert Nullstellensatz(see |Kot Corollary 1.7]), 

for some j ^ 11^=1 polynomial T{xY lies in the polynomial ideal / generated over K by Gj's. It 

follows that some divisor Q{x) of T{xy belongs to the reduced Grobner basis of / (w. r. t. the "lex" 
order). Since this element Q{x) can be constructed via Buchberger's algorithm, its coefficients belong to 
K. In addition, we have Q{f) = and degQ ^ J degT ^ (0^=1 • 

Remark 2.3.3. Note that the existence of such a polynomial Q follows also from Artin's Approximation 
Theorem [A, Theorem 1.7]. 

Corollary 2.3.4. Let 7 be a PBW operad with k binary generating operations. Then the generating 
series of this operad is a solution of an algebraic equation of degree not greater than 4^ . 
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Proof. The algorithm described above impUes the existence of a system of k quadratic and one hnear 
equation on A; + 1 functions. Then we apply Lemma I2.3.2I □ 



2.4. Operads of subexponential growth. Below we present an application of the above theory for 
non-symmetric operads with a small growth. Recall that a sequence {a„}n^o of nonnegative real numbers 
is said to have a subexponential growth if its growth is strictly less than exponential, that is, for each 
d > 1 there exists C > such that < Cd" for all n > 0. 

Corollary 2.4.1. Let 7 be a non- symmetric operad with a finite Grobner basis of relations. Suppose 
that the growth of dimensions IP(n) is subexponential. Then the generating series Gy = X^n^i dim7{n)z'^ 
is rational. In particular, the sequence of dimensions dim'J'{n) should have polynomial growth [n'^] for 
some integer d. 

The proof of Corollary 12.4.11 is based on the general facts about the Taylor coefficients of algebraic 
functions and on the positivity of coefficients in the system (|2.2.2p . 

Consider the system of equations (|2.2.2p obtained while using the first algorithm. Let us remind that 
the series yi = X]n>o Vi^nZ^ is a generating series of the set of monomials in M;,. and the right hand side 
of any equation from this system has strictly positive coefficients. Therefore, for all i the series yi is 
an algebraic function with nonnegative integer coefficients bounded by the dimensions of 7{n). We first 
explain the technical result on algebraic functions with subexponential growth of coefficients and then 
present a proof of Corollary 12.4.11 

Lemma 2.4.2. Suppose that f{z) := Yln^i fnZ^ is an algebraic function, such that the sequence of 
coefficients form a sequence of nonnegative real numbers with subexponential growth. Then there exists a 
rational number s, integer number m and a pair of constants C_, such that for uq sufficiently large 

(a) the upper bound fn < C^n^ is true for all n > uq, 

(b) the lower bound fn > C-u'^ is true for at least one index n in each consecutive collection {N, N + 
1, . . . , + m — 1} of m integer numbers, where N > uq. 

Notice that the exponent s in the above Lemma 12.4.21 is one less then the so-called Gelfand-Kirillov 
dimension 

Dim[/] := m ^^(^^=0^^) 



Inn 

(see [KLj for more details about GK-dimensions of algebras). Lemma [2. 4. 31 below implies that the integer 
number m is bounded from above by the number of singular points on the unit circle of the function /. 

Proof. Let us remind the following well known theorem about asymptotic of coefficients of algebraic 
function based on the Puiseux expansion near critical points. 

Theorem (see jFl Theorem D]). If f{z) is an algebraic function over Q that is analytic at the origin, 
then its n-th Taylor coefficient fn is asymptotically equivalent to a term of the form 



n 



fn = |3^n'Y.'^^^^+OW'-' 
i=0 

where s is a rational number, t < s, (3 is a positive algebraic number and uoi are algebraic with \uJi\ = 1. 

It follows from the proof in [F, Theorem D] that /3 is an inverse of the radius of convergence. Therefore 
in our case (where the sequence of nonnegative integer coefficients fn has subexponential growth) j3 is 
equal to 1. The numbers oji are equal to the singular points on the unit circle of function f{z). The 
upper bound (a) of Lemma 12.4.21 follows from the upper bound: 

m m 
i=l i=l 

The following simple sub-lemma implies the lower bound {b) in Lemma 12.4.21 and finishes its proof. □ 

Lemma 2.4.3. For any given collection of different complex numbers uji, . . . ,ujm with |a;j| = 1 and a 
given collection of constants Ci, . . . , Gm there exists a constant C such that for all n there exists an index 
k = k{n) G {n, n+l,...,n + m— 1} such that \ Yl^i ^i^il > C*. 
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Proof. The proof is by induction on m (the number of summands). The induction base is trivial since 
the absolute value ICiwf | = |Ci||a;i|'^ = \Ci\ does not depend on n. 

Induction step. Consider an integer number n. The induction hypothesis implies the existence of a 
constant C which does not depend on n and an integer k = k{n) G {n, n + l,...,n + m— 1} such that 

ni// \\/ \ k 



(2.4.4) 



i=l 



ijJi 



> 2C. 



Therefore the absolute value of either the number | Cm+i + YlT^i 
of the sequence {Cm+i + Ta=i Ci (zj^Tt) ) 



fc+i 



or the next element 



should be greater than C since their difference coincide 
with the left hand side of the inequality ()2.4.4p . The obvious equality of the absolute values 

m+l 



i=l 



m+l 



i=l 



n 



finishes the proof of the induction step. 

Lemma 2.4.5. Let f{z) and g[z) he a pair of algebraic functions such that their Taylor expansion at the 
origin has nonnegative coefficients with subexponential growth. Then the Gelfand-Kirillov dimension of 
the product f{z)g{z) is a sum of the Gelfand-Kirillov dimension of multiples: 

D\m[f{z)g{z)] = Dim[/(z)] + D\m[g{z)] 

Proof. Let us denote the n-th Taylor coefficient of the product f{z)g{z) by {fg)n- We have {fg)n = 
J2i+j=n fidj- "^^^ upper bound on Gelfand-Kirillov dimension is obvious: 



Y,{f9)^ < E C+(/)^°'"t^l"'C'+(5)j°""[^l-^ ^ 



i=l 



where C+(/),C+(g) are the constants from the upper bound in Lemma 12.4.21 for functions / and g 
respectively. We use a similar sequence of inequalities valid for all sufficiently large n in order to prove 
the lower bound: 



= fi9j ^ Yl ^ 



i=l 



n 



n 



Dim[/]-l 



n 



Dim[c?]-1 



where {C-.{f),mf) and {C-{g),mg) are the constants and integer numbers from the lower bound in 
Lemma 12.4.21 chosen for algebraic functions / and g respectively. □ 

Finally, we can prove Corollarv l2.4.1[ 

Proof. Let IP be a non-symmetric operad with a finite Grobner basis such that the sequence of dimensions 
dimlP(n) has a subexponential growth. 

Consider a system of equations (I2.2.2p obtained while using the first algorithm. We know that the 
functions yi{z) are algebraic with subexponential growth of coefficients. Lemma [2.4.21 implies the existence 
of finite nonzero GK-dimension of any infinite series y,. Let us reorder the set of unknowns according 
to the value of their GK-dimension. I.e. we suppose that are polynomials that have zero 

GK-dimension; yiQ+i, . . . ,yi^ has GK-dimension ai; ... ; yi^._-^^i, ■ ■ ■ ,yi^ has GK-dimension Or, where 
< ai < . . . < Ur and Ir = N. For any given s and for any given i £ {h-i + 1, • • • , ^s} the i'th equation 
in System [2.2.21 should be of the following form: 



(2.4.6) 



yi= Yj Pji^)yj + M^^y^^ 

j = ls-l + l 
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,yis 



where Pj{z) are polynomials and fi is a polynomial in (Ig-i + 1) variables with non-negative coefficients. 
Namely, fi depends only on the first Ig-i unknown variables and do not depend on y^'s with i > Ig-i- 
Indeed, Lemma 12.4.51 implies that the GK-dimension of the product of functions yi^ ■ ■ ■ yi^ should be 
greater than as if there exists at least one ij > Ig or if there exists at least two different multiples yi. and 
yi., with ij > Is-i and ij/ > Iq. 

It remains to show that the solutions of the system ()2.4.6p is rational. Actually, it is well known in 
the theory of generating functions that the solutions of a linear system of equations with polynomial 
coefficients are rational functions (see, e.g., [SI]). However, our system is not generally linear, so, we 
need an induction argument. 

By the induction on s, for the vector = . . . ,yij'^ we get the system of the form 

Ys = AsYs + Bs, 

where Ag G Mat(;^_i^_-^)x(is-is-i)(^^[-^]) ^^'^ is a vector of rational functions which are equal to zero 
at the origin. The vector Bg is obtained by substitution of solutions ?/i's with i < Ig-i which should be 
rational by induction hypothesis. Then 

= (Id -As)-^Bs, 

so that all infinite series yi are rational functions. Then the function Gj> = ^ j/j is also rational. □ 

3. Symmetric and shuffle operads 

3.1. Generating series for a shuffle composition. The first change we should do in the case of 
symmetric operads (compared to what we have explained for non-symmetric operads) is to change the 
type of generating series. Suppose that a subset Ai C B{T) defines a monomial basis of a shuffle operad 
V := -F(T)/(<I>) (meaning that V is the quotient of the free operad T(T) by the ideal generated by a 
subset ^ C J-"(T)). The exponential generating series of the dimensions of V is defined as follows: 

E^n ^ar(v) 
dimV(n) — - = > — r^, where ar(v) means the arity of v. 
^ ^ n\ ^ ar(v)\ ^ ' 

If there is an additional grading of the set of generators T such that all relations from $ are homogeneous, 
one can also consider an exponential generating series in two variables. Let = H V{n) be a 
homogeneous basis of V{n) and let M.ri,k be the subset of M.n consisting of the elements of degree k. 
Then we define 



One can equivalently define 



dee m 



TV. — ' — ' n! 



where -fr-p(„)(t) is the Hilbert series of the graded vector space V{n). As well as in the non-symmetric 
case (Section 12. we do not need the additional parameter t in the most of our examples. We provide 
our proofs mostly for the one-variable series meaning that they remain valid also for the two-variable 
ones with minimal additions. 

Similar to the case of non-symmetric operads, one can define a shuffle composition of vector spaces 
li{Vi, . . . ,T'm)sh (where fi is an element of a free shuffle operad 3" and Vi, ■ ■ ■ ,Vm are the graded vector 
subspaces of 3~) as the vector space generated by all possible shuffle compositions 

(3.1.1) /^(Pi, • • • ,?'m)s/i, 

where pi belongs to the graded component Vi{ki) for all i = 1, . . . , m with ki + ■ ■ ■ + km = n. Similarly to 
the non-symmetric case, the set of underlying internally labeled trees in ^("Pi, . . . ,'Pm)sh should have fi 
as a root vertex and the i-th subtree should belong to the basis of Vi . The main difference with the non- 
symmetric operads concerns the external labeling. As it was mentioned in the definition of divisibility in 
Section [L3l a possible external labelings of a tree from fiiVi, . . . ,Vm)sh should preserve the local order 
of minima of leaves in subtrees (see the proof of Lemma 13.1.41 below) . 
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Let R = Q[[z]] be the ring of formal power series. Define a multilinear map C : i?" — )• i? as follows: 

(3.1.2) C{f,g){z):= r f'{w)g{w)dw forn = 2 

Jo 

(3.1.3) and C(/i,...,/„) :=C(/i,C(/2,...,/„)) forn > 2. 
The next Lemma establish a connection of this operation with shuffle composition. 

Lemma 3.1.4. Let fi,Vi, . . . ,Vm be as above and let S = /^(T^i, . . . ,Vm)sh- Then 

Esiz) = C{Er„...,ErJ. 

Proof. By linearity, it is sufficient to check the above relation in the case of one-dimensional vector spaces 
Pi, ... , Vm- Assume that Vi is spanned by the basis element pi of arity n^. Let n = ni + • • • + n^- Then 



Z 

Es{z) = — c(ni,...,n„ 



where c = c(ni, . . . , n^) is equal to dimS'(n). For each A; = 1, . . . , m, denote by N}^ the n^-element set 
{ni + • • • + Uk-i + 1, . . . , ni + • • • + rifc}. It follows from the definition (cf. [DK, Def. 2]) that the number 
c(ni,...,nm) is equal to the number of permutations a € Sn such that mincr(A^i) < mS.\ia{N2) < 
••• < rain a (Nm) and the restriction of a to every Nf^ is an isomorphism of ordered sets. Therefore, 
c{ni, . . . ,nm) is equal to the number of decompositions [l..n] = Qi U ■ ■ ■ U Qm with \Qk\ = Uk and 
minQi < • • • < rainQm (here Qk = (jN^ for some a as above). The first inequality is equivalent to the 
condition 1 € Qi, hence for every Qi3 \ (there are {^~_^ ways to choose it) there is exactly c(n2, . . . n^) 
ways to get decompositions (^2 U • • • U Qm of the same kind for the set \ Qi. Thus, we have the 

relations 

c{ni,n2) = {^^"^ ^ ] and c(ni, . . . , n^) = ( M c(n2, . . . , n^)- 
V ni - 1 y \ni-lj 

For the generating functions, we obtain the equalities 

n j— c(ni,n2) = ni — - — -, 



n—c{ni, ... , Um) = ni — - ■ r7c(n2, . . . , n^) 

which are equivalent to the desired integration equalities. □ 

Remark 3.1.5. The equation in Lemma '3.1.4\ is equivalent to the following system of ordinary differential 
equations for the functions hk{z) = £'^(-p^ ^^^^^(z); 

h[{z) = E!p^{z)h2{z), 
h'^iz) = E!p^{z)h3{z), 

with the initial conditions hk[Qi) = for ^ j ^ m — k. This system uniquely determines the functions 
hi, ... , hrn-l- 

The following easy verified property of the operation C will be used later in Theorem I3.3.2[ 
Proposition 3.1.6. One has C{f,g) + C{g,f) = fg and, generally, 

C{fa{l),- ■ ■ , fain)) = /1/2 ■ ■ ■ fn- 

cr&Sn 

In particular, C{f, ...,f) = ^. 

In the view of Lemma 13.1.41 this means that the sum of shuffle compositions of some vector spaces 
with respect to all orderings is equal to their non-symmetric composition of the same arity. 
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3.2. System of differential equations. So far we were not able to formulate any statement about 
generating series of an arbitrary shuffle operad with a finite Grobner basis. To establish some properties 
of these series, we require additional assumptions the main of which is given in Definition l3.2.1l below. This 
assumption hold in a number of examples some of which are discussed below. We have checked also that 
for all known symmetric PBW operads there exists a monomial shuffle operad with the same generating 
series as the initial PBW operad but with the following property being satisfied (see Conjecture 14.1.21 
below) . 

Definition 3.2.1. • The planar skeleton of a shuffle monomial m in a free shuffle operad is the 

corresponding planar internally labeled tree, that is, it is obtained from m by erasing the labels 
(numbers) of all leaves. 

• A subset Ai of monomials in the free shuffle operad ^(T) is called shuffle regular if for each 
monomial m G all monomials with the same shuffle skeleton belongs to M. 

For example, the set 

a(/3(xi, X2), 7(x3, X4)), a(/3(j;i, X3), 7(x2, 2:4)), a(/3(xi, 2:4), 7(^2, X3)) 

forms a shuffle regular subset with a shuffle skeleton a(/3(,),7(,)). 

• A monomial operad V is shuffle regular iff the corresponding monomial basis is a shuffle regular 
subset. 

It is obvious that a monomial operad is shuffle regular iff the set of generating monomial 
relations is shuffle regular. 

• Given a set of generators T of a symmetric or shuffle operad V and an admissible ordering of 
monomials, the operad V is called shuffle regular if the set of leading terms of the corresponding 
monomial ideal of relations is shuffle regular. In other words, there exists a reduced Grobner 
basis of the ideal of relations of V in J^(T) with shuffle regular set of leading terms. 

Example 3.2.2. According to the Grobner bases calculated in [DKj . one can see that the operads Com, 
AntiCom and Assoc are shuffle regular [DK, Examples 8,10] whereas the operads Lie and PreLie are 
not (with respect to given orders on shuffle monomials) [DKl Examples 9,11]. On the other hand, if we 
change the ordering of monomials by the dual path-lexicographical ordering, the both operads Lie and 
PreLie become shuffle regular. 

For instance, consider the shuffle operad of associative algebras Assoc. The leading terms of a Grobner 
basis of the ideal of its relations are listed in \DK\ Example 10]. They are the shuffle monomials with the 
shuffle skeletons 

a(a(ai, -)-), a(/3(ai, -)-), /3(/3(ai, -)-), 

where a : a{x, y) = x-y and /3 : /3(x, y) = y- x are the generator operations for Assoc. The cases of other 
operads listed above are analogous. 

Theorem 3.2.3. LetV be a shuffle regular symmetric operad such that the corresponding set of generators 
and a Grobner basis of relations are finite. Then there exists an integer N and a system of integral 
equation on N + 1 functions yo = yo{z, t), . . . ,yiy = yN{z, t) 

(3.2.4) y, = i'^' QiC{ys„...,ys,J for i = 1, . . . , N, 

se[o..N]'^i 

such that Ep{z, t) = yi{z, t), yo = z and yi{0, t) = yi{z, 0) = for all i > 0. The numbers q\ G {0, 1} 
and the nonnegative integer numbers di, ai and N are bounded from above by some functions of the degrees 
and the numbers of generators and relations of the operad V . 

Our proof of Theorem 13.2.31 (as well as the proof of Theorem 13.3.21 below) is close to that of Theo- 
rem I2.2.1[ The main difference is in the counting the number of external labels of a planar tree. This 
reduces to a simple change in the right-hand side of the formula (j2.2.2p : 

Vsi ■ ■■■■Vsd, ^ C{ys^,...,ys^^). 

Namely, one should replace the product of the series by the sign of the operator C(. . .) applied to them. In 
order to make our exposition in symmetric case more self-contained we repeat one of the proofs-algorithms 
in all details. 
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Proof of Theorem \3.2.3l Suppose that an operad V has a finite set of generators T and a finite set 
of monomial relations (It is enough to consider the monomial relations since we are dealing with 
generating series, therefore there is no difference between the relations that form a Grobner basis and 
the monomial relations presented by the leading terms of the first ones.) Let d be the maximum level 
of leaves of elements in ^ (by the level of a vertex/leaf in a tree we mean the number of vertices in a 
path from the root to this vertex/leaf). As it was mentioned in Proposition 11.3.1] every monomial u in a 
free operad J'{T) generated by T may be identified with a rooted planar tree whose vertices are marked 
by elements of T and leaves are numbered by natural numbers 1, . . . ,ar{v) in such a manner that this 
numbering preserves the ordering of minimums in each internal vertex. Given such a monomial v, by its 
stump b{v) we mean the shuffle skeleton of its maximal monomial left divisor such that the leaves and 
the internal vertices of b{v) have levels strictly less than d. 

Let Stump be the set of all stumps of all nonzero monomials in V. Let be the cardinality of this set. 
The elements 5i, . . . , 5jv of Stump are partially ordered by the following relation: bi < bj iS i ^ j and bi 
is a left divisor of bj as a rooted planar tree. Let be the set of all monomials in (= the monomial 
basis of) the right-sided ideal generated by all possible versions of the internal labeling of a stump 6j, 
and set 

j:b,<bj 

The sets A^^. have empty pairwise intersections. Moreover, the set Ui^i -^bi form a monomial basis of 
the operad V. We have 

N 

1=1 

where yi{z) = E^^^^^j^^{z) is the exponential generating series of the span of the set Mbi- For every 

element (=operation) ^ G T of some arity n, define the numbers . . . , in) for all 1 ^ ii, . . . , z„ ^ A 

as follows: 

f 0, if /i(A?r,---,Mr)s/, = inP 

(3.2.5) jf,{h,---,in) = { ' ~ 

[ j, a the stump b{n{Mb,^,..., Mb, J) = bj. 

(where the shuffle compositions for monomial sets are defined as the union of all compositions of 
type (I3.1.ip ). Note that the sets of the type fi{M.bi^ ' • • • T^b,„)sh have vanishing pairwise intersections. 

Let u be a nonzero monomial in V with the root vertex labeled by fi. Then v G fi{Mbi^ ) • • • ) ^bi„ )sh for 
some fi,ii, . . . ,in, that is, v = ^{vi^^ . . . ,fj„)o- where the monomials Vi- G M.b,, and a shuffle composition 
a are uniquely determined by v. Suppose that the index j = jf_i{ii, . . . ,in) from (I3.2.5P is different from 
zero. Then the shuffle regularity condition of a Grobner basis and the bound on the level of relations and 
stumps implies that v should belong to A^f, . . We comes up with the following disjoint union decomposition 
for all i = 1 ... A 

(3.2.6) Aibj = 1^ ^(-^6i^ ; • • • ;-^fei„)5/i) where fi is a root vertex of each v G -Mbj- 

Then the equation ()3.2.4p corresponds to the generating functions yi{z) = E^^^^^j^ ^{z) (where aj is the 
value of the corresponding grading on the operation fj,). Similarly to the non-symmetric case, it follows 
that A^o consists of the identity operation and all Aib, should contain elements of positive degrees in 
generators. This imply the initial conditions on series yi. □ 

Corollary 3.2.7. Let V be a finitely presented symmetric operad with a finite shuffle regular Grobner 
basis of relations. Then there exists a system of ordinary differential equations 

n 

(3.2.8) y[{z) + ^ q]^iyjy'i = gi{z) for i = 1, . . . ,n 
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where Qji^Q and gi{z) G Q[z] with the initial conditions yi{0) = ■ ■ ■ = yn(0) = 0, whose unique formal 
power series solution {yi{z), . . . ,yn{z)) satisfies the equality 

Ev{z) = yi{z) H hyAf(z) 

for some N ^ n. 

Proof. Let us introduce the functions hi as in Remark 13.1.51 for all combinations of the power series yi 
which appear in the equations in the statement of Theorem 13.2.31 Then the i-th equation given in the 
statement of Theorem 13.2.31 is equivalent to an equation of the form 

n 

P=l sG[l..Af]P 

(which is obtained by the differentiation) with the initial condition yi(0) = 0. After the re-naming 
yN+j = hj and adding the equations of the form 

VN+j = y'kVi 

with the same initial conditions yAr_|_j(0) = (cf. Remark 13. 1.5p . we obtain a system of equations of the 
desired form which is equivalent (up to the ghost variables yN+j = hj) to the initial system. □ 

Remark 3.2.9. The number of the equations in the system of differential equations in some cases can 
he reduced. To do these, one can apply to the shuffle regular monomial operads the same methods as we 
have discussed in Subsections \2.2.^ and \2.2.'3\ for the non-symmetric operads. The first of these methods 
is illustrated in Example \3. 5.6\ below. We leave the detailed description of the algorithms to an interested 
reader. 

Corollary 3.2.10. The exponential generating series Ep(z) of a finitely presented operad V with a 
shuffle regular Grobner basis is differential algebraic over Qo That is, there exist a number n ^ and a 
non-constant polynomial 9 in n + 2 variables such that 

e{z,Ev{z),E'T,{z),...,E^^\z)) = 0. 

Proof. By Artin's Approximation Theorem for differential equations |DH Theorem 2.1], for each positive 
integer a there exists another power series solution {yi{z), . . . ,yn{z)) of the system (j3.2.8p such that all 
functions yi are differentially algebraic and 

Viiz) = yi{z) mod z"- 

for all z = 1, . . . ,n. Taking a = 1 and using the fact that the solution with a zero constant term is 
unique, we conclude that yi{z) = yi{z) for all i. Since the sum of differential algebraic functions is again 
differential algebraic, we conclude that E-p^z) = yi(z) + • • • + yN{z) does satisfy a differential algebraic 
equation. □ 

3.3. Relation sets with additional symmetries. In some cases there exists even more symmetries of 
a given Grobner basis of a symmetric (or shuffle) operad. We present below an additional definition of 
symmetric regular operads. This condition implies even more restrictions on the set of generating series 
(Theorem 13.3.21 and Corollary 13.3.41 below) . 

Definition 3.3.1. • The tree skeleton of a shuffle monomial m in the free shuffle operad J- is the 

corresponding internally labeled rooted (non-planar) tree, that is, we erase the labels (numbers) 
of all leaves and forget the planar realization of the labeled tree m. 
• A subset M of monomials in the free shuffle operad J-{V) is called symmetric regular if for each 
monomial m G all monomials with the same tree skeleton belongs to M. 

For example, the set 

f a(/3(xi, X2), 7(2^3, 2^4)), a(/3(xi, xs), 7(x2, X4)), q(/3(xi, X4), 7(x2, xs)), 1 
\ a(/3(x2, X3), 7(xi, X4)), a(/3(x2, X4), 7(xi, X3)), a(/3(x3, X4), 7(xi, X2)) J 

forms a symmetric regular subset with a tree skeleton a(/3 (-,-), 7(-, -)). 

^If we consider power series on two variables t and z, then the coefficient ring Q is replaced by the ring Q[t] with the 
trivial differentiation -^t = 0. 

az 
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• The definitions of symmetric regular monomial operad and arbitrary symmetric regular operad is 
analogous to the ones given in Definition 13.2.11 for shuffle regular case. 

Obviously, the standard monomial basis of a symmetric regular operad is again symmetric regular. 

Theorem 3.3.2. If the set of leading terms of a finite Grobner basis of a shuffle regular operad V form 
a symmetric regular set then the corresponding system of recursive differential algebraic equations (|3.2.4|) 
reduces to the system of algebraic equations 

(3.3.3) = qi-ys,-...- ys,^ for each i = l,...,N 

* se[0. .m]'^i 

for some formal power series yi, - ■ ■ ,yN with non-negative coefficients such that Ep(z) = miyi + • • • + 
rrinyN for some integers mi, . . . , m„. 

Proof. Consider the algorithm given in the proof of Theorem lS. 2.31 applied to a symmetric regular operad 
v. Note that the generating series of M.b and M-i,/ coincide if the stumps b and b' has the same tree 
skeleton. Moreover, for each collection of monomials pi, . . . ,pn and each n-ary operation fj, and a per- 
mutation a £ Sn there is a bijection between the tree skeletons of the elements of the set fi{pi, ■ ■ ■ ,Pn)sh 
and the ones of the elements of the set fJ.{pcr{i), ■ ■ ■ iPa{n))- Consider a relation from the system (|3.2.4|) 
which corresponds to a given stump 6j: 

s<^[O..NY 

While changing the subtrees of the root operation in a shuffle monomial one may change the planar 
skeleton, whereas the tree skeleton remains the same. Therefore, 

yi = t"- ^(y,^^^),...,!/,^^,))) = — Y QlVsi ■ ■ ■ ■ ■ Vs^ 

se[O..N]'i ' cr&Sd ' s£[O..N]'^ 

The last equality follows from Proposition 13 . 1 .61 Thus, the system (j3.2.4p of integration relations can be 
replaced by the system of algebraic equations. Moreover the equations are numbered by the appropriate 
subset of tree-skeletons. Note that these algebraic equations are numbered by the tree skeletons of the 
monomials whose levels are less than the maximal level of the relations. □ 

The above statement is illustrated in Example 13. 5. 6[ 

Analogously to the case of non-symmetric operads, the classical elimination theory implies the following 

Corollary 3.3.4. The exponential generating series E-p of a symmetric regular finitely presented operad 
V is an algebraic function. 

3.4. Operads of restricted growth. We present here an application of the above theory to symmetric 
and shuffle operads of a restricted growth. We say that a sequence {a„}n^o of nonnegative real numbers 
has a subfactorial growth if for each positive constants A,B > there exists a constant C > such that 
ttn < C{^)b for all n sufficiently large. In other words, this means that the growth [a„] of this sequence 
is less than the growth [nl]. In particular, if the sequence is bounded by an exponent C" then its growth 
is subfactorial. 

Corollary 3.4.1. Let V be a symmetric or shuffle operad with a shuffle regular finite Grobner basis. 
Suppose that the growth of the sequence of the dimensions dimV{n) is subfactorial. Then the exponential 
generating series E-p should satisfy a linear differential equation with constant coefficient. Equivalently, 
the usual generating series G-p = ^„>i dim'P{n)z'"' should be rational. In particular, the sequence of 
dimensions dimV{n) should have an exponential or polynomial growth with integer exponent. 

The proof is similar to the one of Corollary 12.4.11 The key point is to reduce the system of equa- 
tions (|3.2.4p for generating series to a system of linear recursive equations on coefficients. 

First, let us discuss some general conclusions of the subfactorial growth. The series = X]„>o yb,n^ is 

an exponential generating series of the set of monomials in M^, hence the right hand side of any equation 
from the system (|3.2.4|) has nonnegative coefficients. Therefore, for each stump h the coefficient y;, „ of 
the exponential series yi, should be a nonnegative integer bounded by the dimension J'(n). Each equation 
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is numbered by an appropriate stump. However, the proof presented below uses rather the type of the 
system than the combinatorics of stumps and the statement of Coronarv l3.4.1l is still true in a much more 
general case arising in some areas of combinatorics. 

There is a standard combinatorial data attached to the system (j3.2.4p . (See, e.g., \FS\ p. 33]). We say 
that a stump b depends on a stump b' if the right hand side of the corresponding recursive equation for 
generating series yb contains a nonzero summand of the form C{y si, ■ ■ ■ ,ysm)^ where some Sk is equal to 
b'. In other words, the right hand side of the equation (|3.2.4p corresponding to yi, depends on y^' in a 
nontrivial way. We say that the dependence is nonlinear if the right hand side of the recursive equation 
for yh contains a multiple C{ysj, . . . ,ysj., . . . jj/s^), where ys^ = yt and at least one of the series ys^ for 
j ^ k is infinite. 

Let us define a graph of dependencies for a system of recursive equations ()3.2.4p . It is a directed graph 
with vertices numbered by all possible stumps. A pair of stumps b and b' is connected by an arrow if 
the stump b depends on b'. This graph is called the dependence graph and will be denoted by TiV). 
Whereas that dependence graph do not contain all information about the system, it sometimes gives 
growth conditions which we illustrate in the lemma below. 

Lemma 3.4.2. Given an edge b ^ b' in the dependence graph T{V), there exists an integer d and a 
polynomial a{n) with positive integer values for all sufficiently large integer n such that the following 
inequality is satisfied for the coefficients of the corresponding generating series: 

(3.4.3) yb,n ^ a{n)yy^n^d- 

Moreover, if the dependence b ^ b' is nonlinear one may choose a polynomial a{n) different from constant. 

Proof. By definition (see (I3.2.6P for details), there is an arrow 6 — )• 6' in the dependence graph if and only 
if there exists a collection of stumps (6i, . . . , br) such that b' = bj for some j and there is an embedding 
of sets 

IJ.{Mbi,...,Mbr)sh cMb- 

For each i ^ j let us choose an element Vi G J^bi- Let di be the arity of the corresponding monomial 
Vi £ V. Consider the subset 

Mb fe' ■= I fJ,{wi,...,Wj-i,Mb',Wj+i,...,Wr)sh \ c u(Mb Mb)sh 

^ \ s.t. Wi has the same shuffle skeleton as Vi for all i ^ j j ^' ' 

The number of elements of arity n in the set Aib^b' should not be greater than the number of elements 
of the same arity in A4b- The detailed counting of the elements in Aib^b' using Lemma 13.1.41 provides 
the inequality 



/ n 






[di 











X 



dj~i - 1 /V Ei>j di . 



which shows the existence of d and a{n) as prescribed in the lemma. Here each binomial coefficient in the 
left hand side is equal to number of shuffle monomials in the corresponding set Mbi with shuffle skeleton 
Vi. Thus, if the dependence 6 — )• 6' is nonlinear, then there exists a collection of monomials Vi such that 
the corresponding product of binomial coefficients is different from constant. □ 

Lemma 13.4.21 has a simple corollary for operads with small growth of dimensions: 

Corollary 3.4.4. If the growth of dimensions V{n) is subfactorial then any loop 6i 62 &i 
in the dependence graph r('P) does not contain nonlinear dependencies. 

Proof. Consider a collection of pairs [(di,ai(n)), . . . , {di,ai{n))\ such that the sequence of following in- 
equalities holds: 

/ I .-1 \ 

2/f>i,n. ^ ai{n)yb2,n-di ^ ai{n)a2{n - di)yb^,n-d ^ • • • ^ 11 "j^'^ ~ X] yhi,n-(di+...+d„)- 



1=1 
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Suppose that there is a nonhnear dependence in the given loop. Then the degree of the polynomial 
a{n) := ^11^=1 ~ Si=i "^j)) is positive, so that a(n) > Cn*^ for some /c ^ 1,C > and for all 

n sufficiently large . Therefore, ybi,n ^ Cn^yh^^n-d ^ C™n''(n - d)^ . . . {n - {m + l)d)''ybi,n-md for 
each m ^ n/d. In particular, if the series y^-^ is different from zero there should exists an infinite 
arithmetical progression of indices {no, no + d, no + 2d,...} such that the corresponding sequence of 
coefficients {ybi,no^ybi,n+o+d, ■ ■ •} does not contain zeros. Put m = [^J • Then 

where A = 2/C and B = 2d+l. Therefore, the only chance for y^^ to have a subfactorial growth is to have 
all polynomials aj(n) equal to positive constants. In particular, all dependencies should be linear. □ 

Proof of Corollary \3.4.1\ The proof is by induction on the number of possible stumps or by the number 
of vertices in the dependence graph T{V). The induction base easily follows from Corollary 13.4.41 since 
any arrow in a graph with one vertex should be a loop. Therefore, the recursive equation on a coefficients 
of generating series should be linear. 

Induction step. Let F be a maximal proper subset of vertices in a graph T(V) such that there is 
no outgoing arrows to the remaining set of vertices V. Let G be a maximal subgraph spanned by V 
respectively, let G be a subgraph spanned by V and all arrows between them. (We omit arrows coming 
from V to V.) Notice that the subgraph G may be empty and, on the contrary, G should contain at least 
one vertex. Moreover, the subgraph G should be a dependence graph for the subset of stumps numbered 
by vertices in V. From the induction hypothesis it follows that for each b € V the corresponding usual 
generating series Gb = X]„>o yb,nZ^ is rational. Lemma 13.4.61 given below shows that any summand of 
the form C{yb^, . . . ,ybi) where all bi belong to V should be an exponential generating series of a sequence 
such that the corresponding ordinary generating function is rational. On the other hand, the maximality 
property of V implies that each two vertices from V should be connected by a directed path. Therefore, 
any arrow b ^ b' G G should belong to a loop where the remaining part of the loop is a directed path 
from b' to b. Corollary 13.4.41 implies that all dependencies in this wheel should be linear. Hence, the 
system for the usual generating series Gb with b £ V reduces to a system of linear equations with rational 
coefficients. As we have seen in the end of Section [2.41 this implies that all ordinary generating series are 
rational. □ 

The following Lemma is well known. We include its simple proof for completeness. 

Lemma 3.4.5. Let G{z) = X^j>i cin-z" and E[z) = X]i>i ^-^^ ordinary and exponential generating 

functions of the same sequence of complex numbers {a„}„^i. Then the function G{z) is rational if and 
only if the function E{z) satisfies a non-trivial linear differential equation with scalar coefficients. 

Proof. The condition "G(z) is rational" means that a recurrent equation 

k-l 
3=0 

holds for all n ^ 1. It is equivalent to the recurrent relation 

{n + k)\ _^ (n+j)!, 

On+k - l^Cj — bn+j 

j=0 

for the numbers bn = an/nl. This is equivalent to the differential relation 



i=o 

for the generating function E{z) = Yln>i ^n-^"- I— ' 

Let £ be the set of all exponential generating series such that the corresponding ordinary generating 
series are rational functions. In other words, £ is the set of all exponential generating functions which are 
the solutions of non-trivial linear differential equations with scalar coefficients. Now, the next Lemma is 
obvious. 

25 



Lemma 3.4.6. The set £ is closed under multiplication, differentiation and integration. In particular 
if Ef,Eg £ E then C{Ef,Eg) = E'j{w)Eg{w)dw also belongs to £, that is, the corresponding ordinary 
generating series is also rational. 

3.5. Examples for symmetric operads. 
3.5.1. Examples of shuffle regular PBW operads. 

Example 3.5.1. Consider the class of so-called alia algebras introduced by Dzhumadil'daev [Dzh], that 
is, the algebras with one binary operation (multiplication) and satisfying the identity 

{[xi,X2],a;3} + {[x2,a;3],xi} + {[x3,xi],X2} = 0, 

where [xi,a:2] = xiX2 — X2X1 and {xi,X2} = X1X2 + X2X1 (these algebras are also referred in [Dzh] as 
1-alia algebras). 

Let us choose the generators a : (xi,X2) ^ [xi,X2\ and /3 : (xi,X2) ^ {xi,X2} for the corresponding 
symmetric/shuffle operad Alia. Then the leading term of the relation corresponding to the above identity 
(with respect to the path-lex order with /3 > a) is /3(xi, a(x2, X3)), that is, the only shuffle monomial 
corresponding to the shuffle skeleton /3(— , a(— , — )). 

Obviously, there is no overlapping of the leading term /3(xi, a(x2, X3)) of the relation with itself, hence 
the relation is the unique element of the Grobner basis of the relations of Alia. Then we have the 
following three elements of the set of stumps B (in terms of the proof of Theorem I3.2.3p ): 

Bo = Id , = a, B2 = p. 

We get the relations 

Hi = C{Ej^iia,E_Aiia), 

y2 = C{EMia,yo + y2), 

EAlia =2/0 + 2/1+2/2 

(here we use the linearity of the operation C to shorten the summation in the right hand sides). Using 
the linearity of C and Proposition 13.1.61 we get the system 

2/0 = z, 
2/1 = E%J2, 

y2=E%j2-C{EMia.yi), 
EAlia = 2/0 + 2/1 + 2/2, 

which leads to the equation 

C(y,yV2) =z-y + y2 
for y = EAiia{z). After differentiation, we get the equation 

2/'yV2 = 1 - 2/' + 22/2/'- 
Using the initial condition y(0) = 0, we get the algebraic equation 

2/^/6 -y^ + y = z 

for y. In particular, the function y = EAiia{z) is algebraic. 

Since the operad Alia has a quadratic Grobner basis of relations, it follows from [Hj (see also [DKl 
Cor. 3]) that the operad Alia is Koszul. By the Ginzburg-Kapranov relation, its exponential generating 
series y satisfy the relation 

f{-y) = -z, 

where f{z) is the exponential generating series of the quadratic dual operad Alia' . It follows from the 
equation above that 

Ej^ua'\z) = z + z' + z'^/Q. 

This polynomial coincides with the result of a direct calculation given in the concluding remark of [Dzh]. 

Recall now two another classes of algebras from [Dzh] . A nonassociative algebra is called left (respec- 
tively, right) alia if it satisfies the identity 

/(xi,X2,X3) = [X1,X2]X3 + [x2,X3]xi + [x3,Xi]x2 = 0, 
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or, respectively, the identity 

r{xi,X2,X3) = Xi[x2,X3] +X2[x3,Xi] +X3[xi,X2] = 0. 

Consider the left hand side r{xi,X2,X3) of the last identity. Using the above generators a and /3 with the 
substitution 2xiX2 = a(xi, X2) + P{xi, X2), we see that the leading monomial of r{xi,X2,xs) with respect 
to the same path-lex order with ^ > a is the same monomial l3{xi,a{x2,xs)) as for alia algebras. By the 
same reasons as above, we see that the opcrad of right alia algebras is PBW and Koszul with the same 
generating series as the operad Alia. By the right-left symmetry, the same is true for left alia algebras. 
Thus, we get 

Proposition 3.5.2. The three operads for alia algebras, left alia algebras and right alia algebras are 
Koszul with the same exponential generating series y = E-p(z) satisfying the equation 

y^/6 -y^ + y = z, 

that is, 

2 11 s 25 4 127 . 259 . 1475 7 17369 o 943855 □ 2906189 ^ ^/ iin 

v(z) = z + z^ -\ z-^ H z^ -\ z^ H z^ H z^ H z^ H z^ H z^" + Ofz^M. 

y^^j ^^^^g^^g^^l2 9 18 72 1296 1296 ^^^^ ' 

Each of their three Koszul dual operads is finite- dimensional and has exponential generating series 

E-p^{z) = z + z^ + z^/Q. 

3.5.2. Examples of symmetric regular operads. In the next two examples, we consider the operad of upper 
triangular matrices over non-associative commutative rings. 

Example 3.5.3. Let i? be a commutative non-associative ring (or a k-algebra). Then it is easy to see 
that the algebra UT2{R) of upper triangular 2 x 2-matrices over R satisfies the identity 

[a;i,a;2][x3,X4] = 0, 

where [a, b] = ab — ba. 

Let us denote by NU2 the operad generated by the operation of non-symmetric multiplication : 
(xi,X2) !->■ X1X2 (i. e., the arity two component NU2 is spanned by fj, and fj,' : {xi,X2) i->- X2X1) subject 
to this identity. Consider the corresponding shuffle operad J\fU2 generated by two binary generators, 
namely, the operations /x and a : {xi,X2) i-> [a;i,a;2]- Then the above identity is equivalent to the pair of 
shuffle regular monomial identities 

/i = //(a(-,-),a(-,-)) = and /2 = a(a(-, -), a(-, -)) = 0. 

Therefore, the ideal of relations of the shuffle operad AfU2 is generated by the following six shuffle 
monomials obtained from /i and /2 by substituting all shuffle compositions of four variables (which we 
denote for simplicity by 1,2,3,4): 

mi = ^(a(l,2),a(3,4)), m2 = fJ.{a{l,3),a{2,4)), ms = ^(a(l, 4), a(2, 3)), 
m,4 = a{a(l, 2), a(3, 4)), 77x5 = a(a{l, 3),a{2, 4)), mg = a(a(l, 4), a{2, 3)). 

Let us describe the set B of all stumps of all nonzero monomials in AfU2- Since the relations have their 
leaves at level 2, B includes all monomials of level at most one, that is, the monomials 

So = Id, Si =^(-,-),fi2 = «(-,-). 

For the corresponding generating series t/j = yi{z) with z = 0, 1, 2 we have 

yo = z, 

< yi = C{yo,yo) + C{yi,z) + C{z,yi) + C{y2,z) + C{z,y2) + C{yi,yi) +C{yi,y2) + C{y2,yi), 
, y2 = C{yo, yo) + C{yi,z) + C{z, yi) + C{y2, z) + C{z, 2/2) + C(yi, Z/i) + C(yi, 2/2) + C(y2, yi). 

We see that yi[z) = y2{z) and Ej^i(,^[z) = y[z) = yo{z) + yi{z) + 1/2(2;) = z + 2yi{z). The second equation 
of the above system gives, after differentiation, the ordinary differential equation (ODE) 

y[ = z + 2zyi + 2zy[ + 3yiy[ 

with the initial condition yi{0) = 0, which is equivalent to the ODE 

{y'{z)-l){2-z-3y{z))=4y{z) 
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on y{z), again with the initial condition y{0) = 0. It follows that 
Emu2{z) = y{z) = \{2-z- 2Vl-4z + z2' 



= ^ + ^2 + 2 z3 + + f ^5 + 2|1 ^6 + 413^7 + 2(Km ^8 + 3^9 ^9 + 3^^10 + ^(^10) 

Let us generalize Example l3.5.3l to the case of matrices of order n. The following description of identities 
easily follows from the fact that the diagonal elements of the commutator of two upper triangular matrices 
are zero. 

Lemma 3.5.4. Let Un{R) be the algebra of upper triangular matrices of order n over a (non-associative) 
commutative ring R. Then for each n-ary multiple composition f of multiplications of matrices, the 
identity of 2n arguments 

f{[xi,X2], . . . , [Xn^l,Xn]) = 0, 

holds in Un{R), where, as usually, [a, b] stands for ab — ba. 

For example, for n = 2 we get the only identity [xi, X2] [xs, X4] as above. For n = 3 we have the identities 
fi{[xi,X2], [x3,X4], [x5,a;6]) with i = 1,2, where /i(ai,a2,a3) = (0102)03 and /2(oi, 02,03) = 01(0203) (ah 
other identities are obtained from these two by permutations of variables). 

Consider the operad NlAn of upper triangular matrices generated by the non-symmetric operation of 
multiplication ^ : (xi,X2) ^ X1X2 subject to all these identities. Consider a natural generators of the 
corresponding shuffle operad a : {xi,X2) ^ [xi,X2] and /3 : (xi,X2) ^ X1X2 + X2X1. We immediately see 
that the set of shuffle relations of this operad is spanned by monomials and is symmetric regular. 

Corollary 3.5.5. For all n ^ 2 the exponential generating series Ej^y^ is algebraic. 

Example 3.5.6. Let us find a relation for the exponential generating series for the operad V = MlA-^. 
To do this, we use the appropriate version of the method used in Subsection 12.2.21 

Subject to the binary generators a, (3 the minimal sets of the monomial relations of this operad consists 
of the monomials with one of the following 4 tree skeletons: 

C(C(a(-, -),«(-, -)),"(-,-)), where C,C G 

Then the set T(V) from Lemma l2 . 2 . 71 consists of the following five tree skeletons: 

/ (the identical operation), o := q(-, -), 5 := /3(-, -),A := a{a{-, -), a(-, -)), B := /3(a(-, -), a(-, -)). 

We get the following equations for the corresponding generating functions ya, yb-,yA-,yB-, and yi = E-p: 

yi = z + ya + yb, ( yi = z + 2ya, 

ya = yb = \y]-\{yAya + yBya + yayA + yays) + \y\, <^ < ya = \y'j -2yayA + \y\, 
yA = yB = \yl - \ {yAya + yayA) + ^y\- [ = ^yl - yayA + \y\- 

After elimination of the variables and we obtain the following algebraic equation for yj = Ej^jj.^: 
7//^ + (12 z - 24) y/ V (30 + 8 z + 8O) y/ V (-36 z^ + 24 - 32 z - 64) y/ +9 / - 8 + 16 z ^64 z = 0. 
It follows that 

, o 4 . 167 7 26745 o 44045 „ 36969 ^ ^/ ^^^ 

Ej^Us = z + z^ + 2 z^ + 5 z^ + 14 z^ + —z^ + 130 z^ H —z^ + z^ + z^° + O {z^^) . 

4 64 32 8 

3.5.3. Examples of computations via homology. 

Example 3.5.7. Let us define a class of operads where all relations are relations on the commutators. 

Namely, any given finite set of operations T and finite set of linear independent elements ^ C -7^(T) 
defines an operad 

n ■ X-^Tll/r nw Mkl,2;2],X3] + [[X2,X3],X1] + [X3,X1],X2] = \ 

Q^.= T{TU {[-,-]})/ ^ 5(h-],..., [-,-]) for all^G^ J 

In other words, the operad 2$ is generated by the union of the set T and a Lie bracket [-,-]. The first 
identity is the Jacobi identity for the Lie bracket. The set of additional identities in the lower line consists 
of substitutions of the Lie bracket into the arguments of the operations in (I. e., we put the bracket 
[-,-] in each leaf of g). 

Proposition 3.5.8. The above set of relations of the operad Q<j, form a shuffle regular Grobner basis 
with respect to the dual path-lexicographical ordering. (The ordering of generators does not matter.) 
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Proof. There are two types of nontrivial intersections of leading terms of relations in Q$: 

• The leading term [xi, [x2, xs]] of the Jacobi identity intersects with itself. The corresponding 
s-polynomial written for the monomial [xi, [x2, [x3,X4]]] reduces to zero because the operad Lie 
is PBW according to the dual path-lexicographical ordering. 

• For any given 5 G $ the leading term ^([, ],...,[, ]) of the relation ],...,[, ]) of arity m has m 
different intersections with the leading term of the Jacobi identity. 

The typical s-polynomial of the second kind looks as follows: 

S := g{..., [xi, [xj,Xk\], ...) - g{...,J {xi,Xj,Xk), ■ ■ ■), 

where J{xi,Xj,Xk) ■= [xi, [xj,Xk]] — [[xi,Xj],Xk] + [[xi, Xk], Xj] is a notation for Jacobiator. One should 
make the following reductions: 

S + g{..., [[xi,Xj],Xk],.. .) - g{..., [[xi,Xk],Xj], . . ■) - {g - g){. ■ ■ , J{xi,Xj,Xk), . . .) = 

in order to reduce the above s-polynomial to zero. □ 

Starting from a Grobner basis one can define a resolution (also called a DG model) for the monomial 
replacement as it is prescribed in [DKl] and then deform the differential to get a resolution for Q$. In 
order to compute the generating series of the operad it is enough to describe the generators in the last 
resolution. This is given in the following 

Proposition 3.5.9. There exists a resolution of Q$ generated by the union of the following sets: 

• the set of generators T with homological degree zero, 

• the set of right-normalized commutators L := {[xi, [x2^ ■ ■ ■ , [xk^i,Xk] ■ ■ ^ 2}, (The homolog- 
ical degree of the commutator on k-letters is equal to (k — 2).) 

• the union of sets {g{L, . . . , L)sh} for all g £ ^ (The homological degree of g{li, ■ ■ ■ , lk)sh is equal 
to^deg(/i) + l.; 

The generating series of the right normalized commutators L is {e~^ -\-z—\). Multiplying the generating 
series of each of the above sets by (—1)^, where h is the homological degree, we get the equation 

(3.5.10) E^^{z) = 1 - e"^ - Er{z) - ^$(6"^ + z-l) 

for the functional inverse of the exponential generating series of the operad Q$. 

Example 3.5.11. Consider a particular example of the above construction consisting of the class of 
Lie-admissible non-associative algebras. Let V be an operad generated by one binary non-symmetric 
operation (multiplication) subject to the following identities: 

[a, [6,c]] + [6, [c,o]] + [c, [0,6]] =0, 

[a,b][c,d] + [c,d][a,h] = 

where, as usually, [a, b] = ab — ba. The operad V satisfies the conditions of Proposition 13.5.81 with 
T = $ = {/? : (xi,X2) ^ X1X2 + X2X1}. It follows that Eriz) = E^{z) = z'^/2. The identity (13.5. 1U|) 
gives the following equation on the generating series E-p{z) after appropriate substitutions: 

E^\z) = 1 - - ^ - ^ ^ ~ ^3-2Ev + 2Er e"^^ + e'^^^ - 4 e"^^ = 2z. 
It follows that 

, 11 o 49 4 1219 . 811 R 75919 7 97175 o 25827439 „ 116679221 ^ ^, in 

Eviz) = z+z^^ z^+—z^^ z^+ z^+ z^+ z^+ z^H z^°+0 (z^^) . 

' ^ ^ 6 12 120 30 1008 448 40320 60480 ^ ^ 

4. Remarks and conjectures 

4.1. Operads with quadratic Grobner bases. Let us formulate a couple of conjectures about PBW 
operads (i.e. the operads which admits a quadratic Grobner basis). 

Conjecture 4.1.1. The exponential generating series of a symmetric PBW operad with binary generators 
is differential algebraic. 

This conjecture is motivated by the following stronger conjecture that has been checked in all known 
examples. 
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Conjecture 4.1.2. Let V be a symmetric operad generated by binary operations and with quadratic 
relations that form a Grobner basis according to one of the admissible orderings of shuffle monomi- 
als. Then there exists a monomial shuffle regular operad Q with the same dimensions of operations 
(dimV{n) = dimQ[n) for all n). 

Why this conjecture might be true? The most important assumption is that the operad V is symmetric. 
Therefore, one should use the representation theory of the permutation groups $]„,. First, using the 
upper-triangular change of basis one can choose the set of generators T of an operad V such that the 
transposition element (12) G S2 preserves this set. This means that Va G T (12) • a = {—lya' where 
a' G T. Suppose, for simplicity, that (12) - a = ita for all a £ T. For any given pair of generators a, a' the 
set of all operations of the form a o a' form a 3-dimensional representation of S3 which is the induction 
from a one-dimensional S2-representation. Therefore, it decomposes into the sum of one-dimensional and 
2-dimensional irreducible representation. Relations for the 2-dimensional representation looks as follows 



On the other hand, there are two shuffle regular sets of quadratic monomials: 

{a{a'{xi,X2),X3),a{a'{xi,xs),X2)} and {a(xi, a'(a;2, X3))} 

The hope is that the monomial operad Q can be constructed by taking the leading term {a(xi , a' {x2, X3))} 
of the relations that nontrivially projects on corresponding one-dimensional Ss-irreducible subspace and 
the set {a(a'(a;i, X2), 2:3), a(a'(xi, X3), X2)} for those relations which projects on the two-dimensional 
E3-irreducible subspace. 

4.2. Varieties of associative algebras. The symmetric quotient operads of the operad Assoc (i.e., 
the operads corresponding to the varieties of associative algebras with polynomial identities; called here 
associative operads) form, probably, the widest extensively classified and studied class of operads. Note 
that the traditional terminology for these operads differs from the one we use. Particularly, the dimensions 
diuiVn of an associative operad V are usually referred to as codimensions of the corresponding varieties 
of algebras with polynomial identities; the exponential and the usual generating series of the operad are 
referred as the codimension series and the exponential codimension series of the corresponding variety. 



We do not know whether each symmetric quotient operad of Assoc has a finite Grobner basis. On 
the other hand, one can construct examples of shuffle monomial quotients of Assoc which are infinitely 
presented, that is, which have no finite Grobner basis, at least with the standard choice of generators. 
However, the following properties of generating series of associative operads are close to the properties 
of our operads with finite shuffle regular Grobner bases. 

Firstly, in many cases the dimensions diuiVn for such an operad V can be calculated explicitly (see [U] 
IGZl IDrt IBD] ). In such cases, these dimensions are presented as linear combinations of integer exponents 
of n with polynomial coefficients. In particular, the corresponding exponential generating functions are 
differential algebraic, that is, they satisfy the conclusion of our Corollary 13.2.101 

Secondly, the only non-symmetric associative operads are the finite-dimensional ones defined by the 
identity X1X2 ■ ■ ■ Xn = (operads of nilpotent associative algebras); obviously, the generating series of 
such operads do satisfy the conclusions of our Theorem 12.3.11 and Corollary 12.4.11 Thirdly, the sequence 
of dimensions of an associative operad has either polynomial or exponential growth (by the results of 
Regev and Kemer, see [GZl Theorems 4.2.4 and 7.2.2]). This fact is similar to our Corollary 13.4. 1[ 

After these observations, we think that our Corollary 13.2.101 (which suggests Expectation [T]) holds in 
a more general situation, that is, for a wider class of symmetric operads containing both operads with a 
finite shuffle regular Grobner basis and associative operads. 

Conjecture 4.2.1. The exponential generating series of any associative operad is differential algebraic. 

Note that the two languages of universal algebra (based on varieties vs the operads) are equivalent 
provided that the basic field k has zero characteristic. In this case, our Conjecture 14.2.1] is equivalent to 
the following: the exponential codimension series of each variety of PI algebras is differential algebraic. 
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a(a'(xi, X2), X3) -I- a(a'(x2, X3), xi) -|- a{a'{x3, xi), X2) 
Relations for the one-dimensional representation looks as follows: 




(See, e.g., |GZ].) 



A stronger version of Conjecture 14.2.11 holds for associative operads of polynomial growth, that is, 
their ordinary generating series are rational [BDl Ex. 13]. However, for some natural associative operads 
(e.g., for the operad defined by the identities of 3 x 3 matrices) the ordinary generating series are not 
rational [BR]. This means that the conclusion of our Corollary 13.4.11 does not hold for general associative 
operad. In particular, this means that some associative operads have no finite shuffle regular Grobner basis 
(whereas they are always finitely presented, by the Kemer's famous solution of the Specht problem [Kg]). 
This means also that some associative operas are not generic in the sense of our Expectation [3j 

4.3. Detailed description of generating series. It would be interesting to find a closer class of formal 
power series than the class of differential algebraic ones which include all (exponential) generating series 
of generic symmetric operads. For the operads with a finite shuffle regular Grobner basis, we propose 
the following approach. 

Let us define a binary operation on Q[[z]] hy f * g := C{f,g), where C{f,g){z) = f {w)g{w) dw 
is defined in (j3.1.2p . Recall that the multiple operation C from equation (j3.1.3p is also defined via this 
binary operation: C(/i, . . . , /„) = (/i * (/2 * • • • (/n-i * fn) ■ ■ ■ )■ The algebra 'Q[[z\\ with respect to the 
operation f * g satisfy the identity 

{a*h)*c = a*{b*c + c*b), 

that is, this is a Zinbiel algebra (also known as dual Leibniz algebra, or chronological algebrajl. Then 
one can consider the main system (I3.2.4P as a set of algebraic equations over the Zinbiel subalgebra 

QN cQ[[z]]. 

In this connection, the following questions arise. 

1. Does it follows that the exponential series function E'p{z) satisfy a single algebraic equation over 
Q[z] with respect to the operation that is, does it follows that E'p{z) is Zinbiel algebraic over (Q[2:], *)? 

2. What are the conditions for the (non- negative rational) coefficients for a formal power series to be 
Zinbiel algebraic (or for a collection 2/0(2), . . . ,yM{z) of such formal power series to satisfy a system of 
Zinbiel algebraic equations)? 

The answer to the second question could give a detailed description of generating series of, at least, 
operads with shuffle regular Grobner basis. 

References 

[A] M. Artin, Algebraic approximation of structures over complete local rings, Inst. Hautes Etudes Sci. Publ. Math., 36 
(1969), pp. 23-58 

[BD] S. Boumova, V. Drensky, Algebraic Properties of Codimension Series of PI-Aloebras. larXiv:llll. 03191 /1 (2011) 
[BR] W. Beckner, A. Regev, Asymptottcs and algebraicity of some generating functions, Adv. Math., 65 (1987), pp. 1-15 
[ChSch] N.Chomsky, M. P. Schiitzenberger, The Algebraic Theory of Context-Free Languages , Studies in Logic and the 

Foundations of Mathematics, 35, 1963, pp. 118-161 Computer Programming and Formal Systems 
[Dr] V. Drensky, Combinatorial aspects in Pl-rings, Polynomial identity rings, Birkhauser, Basel-Boston-Berlin, 2004, 

pp. 1-130 

[Dzh] A. S. Dzhumadil'daev, Algebras with skew-symmetric identities of degree 3 (Russian), Sovrem. Mat. Prilozh., 60, 

Algebra (2008), pp. 13-31; translation in J. Math. Sci. (N. Y.), 161 (2009), 1, pp. 11-30 
[DK] V. Dotsenko and A. Khoroshkin, Grobner bases for operads, Duke Math. J. Volume 153, Number 2 (2010), 363-396 
[DKl] V. Dotsenko and A. Khoroshkin, Free resolutions via Grobner bases , math, arXiv:0912. 4895 

[DL] J. Denef and L. Lipshitz, Power series solutions of algebraic differential equations, Mathematische Annalen, 267 
(1984), 2, pp. 213-238 

[F] Ph. Flajolet, Analytic models and ambiguity of context-free languages. Twelfth international colloquium on automata, 
languages and programming (Nafplion, 1985). Theoret. Comput. Sci. 49 (1987), no. 2-3, pp. 283-309 

[FS] Ph. Flajolet, R. Sedgewick, Analytic combinatorics. Cambridge University Press, Cambridge, 2009. xiv+810 pp. 

[CZ] A. Giambruno, M. Zaicev, Polynomial identities and asymptotic methods. Mathematical Surveys and Monographs, 
122, AMS, Providence, RI, 2005 

[H] E. Hoffbeck, A Pomcare-Birkhoff-Witt criterion for Koszul operads, Manuscripta Math., 131 (2010), 1-2, pp. 87- 
110 

[Ke] A. Kemer, Ideals of Identities of Associative Algebras, Amer. Math. Soc, Translations of Monographs 87 (1991) 
[Ko] J. KoUar, Sharp Effective Nullstellensatz, Journal AMS, 4, 1 (1998), pp. 963-975 

[KR] D. Krakowski, A. Regev, The polynomial identities of the Crassmann algebra. Trans. Amer. Math. Soc, 181 (1973), 
pp. 429-438 

[KL] G. Krause, T. Lenagan, Growth of algebras and Gelfand-Kirillov dimension. Revised edition. Graduate Studies in 
Mathematics, 22. American Mathematical Society, Providence, RI, 2000. x-|-212 pp. 



'We are grateful to Askar Dzhumadil'daev for pointing us this fact. 

31 



" ~ p 

[LV] J.-L. Loday, B. Vallette Algebraic Operads, xviii+512 pages available at http://math.unice.fr/branov/Operads.html 
[MSS] M. Markl, S. Shnider, J .Stasheff, Operads in algebra, topology and physics, Mathematical Surveys and Monographs, 

96, AMS, Providence, RI, 2002 
[PP] Polishchuk, A., Positselski, L., Qudratic algebras, University Lecture Series, 37, American Mathematical Society, 

Providence, RI, 2005 

[R] E. S. Rowland, Pattern avoidance in binary trees, J. Combin. Theory Ser. A, 117 (2010), 6, 741-758, 
[S] R.Stanley, Differentiably finite power series. European J. Combin. 1 (1980), no. 2, 175-188 

[SI] R. Stanley, Enumerative combinatorics. Vol. 2. Cambridge Studies in Advanced Mathematics, 62. Cambridge Uni- 
versity Press, Cambridge, 1999. xii+581 pp. 

[U] V. A. Ufnarovsky, Combinatorial and asymptotical methods in algebra, Sovr. probl. mat.. Fund, napr., 57 (1990), 
p. 5-177 [Russian] Engl, transl.: Algebra VI, Encycl. Math. Sci., Springer, Berlin 1995, p. 1-196 

[Z] G. W. Zinbiel, Encyclopedia of types of algebras 2010, arXiv:1101.0267y l (2011) 

408 Simons Center for Geometry and Physics, Stony Brook, NY 11794-3636, USA and ITEP, Bolshaya 
Cheremushkinskaya 25, 117259, Moscow, Russia 

E-mail address: anton.khoroshkinOscgp . stonybrook. edu 

Department of High Mathematics for Economics, Myasnitskaya str. 20, National Research University 
Higher School of Economics, Moscow 101990, Russia 
E-mail address: piont@mccme.ru 



32 



